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AB3IIACI 

An  ei|«aipiee  optinization  approach  is  used  to  incorpo¬ 
rate  optiaal  control  into  the  structural  design  process  for 
large  flexible  space  structures.  The  equations  of  notion  for 
an  uncontrolled  systea  are  developed  by  deriving  the  kinetic 
and  potential  energy  for  the  systea  and  then  using  assuaed 
nodes  to  discretize  the  energies.  These  expressions  are  then 
linearized,  the  Lagrangian  foraed,  and  Lagrange's  equations 
written  for  the  systea.  An  existing  optiaal  control  law  is 
incorporated  to  fora  the  equations  of  notion  for  the  control¬ 
led  systea.  A  paraaeter  optinization  technique  is  used  to 
ainiaize  the  aass  of  the  Draper/RPL  Configuration  nodel 
involving  eigenspace  optinization.  A  conputer  algorithm  is 
developed  that  effectively  optinizes  a  global  structural 
paraaeter  vector  to  ainiaize  the  aass  of  the  nodel.  while 
constraining  specified  eigenvalues.  The  eigenvalue 
sensitivities  are  passed  to  a  constrained  function  ainiaiza- 
tion  prograa  called  CONNIN  which  ainiaizes  the  aass  of  the 
appendages.  The  constraints  iaposed  restrict  the  first 
eigenvalue  to  the  left  half  plane  and  the  natural  frequency 
of  the  third  eigenvalue  to  a  specified  stable  region.  The 
result  is  an  algoritha  that  incorporates  an  existing  optiaal 
control  law  into  the  structural  optinization  process. 
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INCORPORATING  CONTROL  INTO  TIB  OPTIMAL  STRUCTURAL 


DESIGN  OF  LARGB  FLBXIBLE  SPACE  STRUCTURES 

I.  imPPB.CII.Qfl 

President  Reagan  has  directed  NASA  to  begin  work  on  a 
permanent  nanned  space  station.  The  Russians  just  completed 
a  record  stay  in  space  and  the  space  shuttle  will  begin 
Bonthly  flights  very  soon.  The  deaand  for  large,  controlled 
flexible  space  strnctnres  will  soon  increase  greatly.  These 
structures  asst  not  only  be  stable,  bnt  aust  also  be  designed 
to  ainiaize  their  aass.  Even  with  the  advent  of  the  space 
shuttle,  transporting  aassive  objects  into  Earth  orbit  is 
costly.  Efforts  to  ainiaize  the  aass  of  the  large  space 
structure  usually  result  in  a  highly  flexible  structure  that 
can  pose  a  difficult  control  problea. 

Current  aethods  of  structural  design  generally  do  not 
incorporate  control  into  the  optiaization  process.  The 
result  is  a  structure  that  is  optiaally  designed  with  respect 
to  strength  or  vibration  or  soae  other  perforaance  paraaeter 
and  then  an  optiaua  control  systea  is  designed  for  the 
existing  structure.  Therefore,  optiaal  design  froa  a 
structural  standpoint  can  represent  varying  structural 
paraaeters  to  ainiaize  the  aass,  whereas  for  a  control 
designer,  optiaization  aeans  controlling  the  nuaber  and  size 
of  controllers  suoh  that  a  specified  perforaance  index  is 
ainiaized  [9]. 


While  many  recent  papers  have  addressed  the  problem  of 
optimizing  a  control  system  for  a  given  space  structure, 
this  paper  will  examine  methods  of  designing  a  large  flexible 
space  structure  with  a  given  control  system  in  mind. 

A  very  stiff  ,  and  hence  a  more  massive,  large  space 
structure  will  require  less  control  energy  than  a  like 
structure  that  has  mueh  less  mass  and  is  more  flexible.  Also 
the  greater  it's  mass,  the  more  expensive  it  is  to  transport 
into  orbit.  Therefore.  there  is  a  tradeoff  between 


minimization  of  mass  of  a  structure  and  the  ability  to 
control  it.  The  minimum  mass  for  a  structure  is  described  by 
Hale  in  [1]  to  be  the  essential  mass  to  accomplish  the 
mission  for  which  the  structure  is  designed.  He  goes  on  to 
say  that  this  structure  can  be  a  collection  of  discrete 
masses  or  be  spatially  distributed.  When  using  discrete 
actuators  on  a  spatially  distributed  system  the  amount  of 
control  energy  expended  can  be  reduced  by  stiffening  the 
structure.  This  effectively  distributes  the  control  forces 
over  the  'essential  structure',  but  it  is  at  the  expense  of 
increasing  the  total  mass  of  the  structure. [1]  Any  variation 
in  the  parameters  must  retain  control abil ity  and  stability 


of  the  entire  structure. 


HISTORY  Of  RELATED  EFFORTS 

Hale  and  Lisowski  in  [1]  develop  an  analytical 
approach  to  the  optimal  structural  design/ control  design 
problem.  They  use  linear  qua  dr  a t i c- c o s t  control  theory  to 
design  an  optimal  control  for  known  structural  parameters.  A 
weighted  total  mass  term  is  added  to  the  cost  functional  and 
the  structural  parameters  are  varied  to  find  the  minimum 
total  cost.  Discrete  equations  are  the  result  of  spatially 
discretizing  the  partial  differential  equations  that 
represent  the  complex  structure  composed  of  many  distributed 
components.  The  calculus  of  variations  is  used  to  derive  the 
necessary  conditions  for  an  optimum  of  the  problem.  The 
result  is  a  hybrid  system  of  coupled  nonlinear  equations  with 
'no  hope  for  an  analytical  solution*.  A  numerical  solution 
composed  of  two  parts  is  discussed.  The  two  parts  »:e:  1) 
find  a  control  design  that  is  optimal  for  given  structural 
parameters,  and  2)  a  procedure  to  update  the  structural 
parameters.  Thus  the  procedure  is  an  iterative  one  and 
because  the  cost  functional  is  linear  in  the  structural 
parameters,  the  second  partial  derivatives  of  the  functional 
with  respect  to  the  structural  parameters  are  zero.  As 
stated  in  the  report,  when  using  this  approach,  care  must  be 
exercised  as  the  extremal  solution  may  only  be  a  relative 
minima  or  maxima.  One  important  aspect  of  this  work  is  that 
nonlinear  effects  are  not  included. 

In  [2]  Hale  and  Lisowski  use  the  same  approach  as  in 


to 


3 


[1]  for  a  reduced  order  problem.  The  transformation  is 
accomplished  using  the  force-free  eigenvectors  to  form  a  new 
basis  for  the  modal  coordinates.  A  new  set  of  equations  are 
formed  to  approximate  the  solution  to  the  original  problem. 
How  well  this  approximate  solution  represents  the  actual 
solution  depends  on  the  participation  of  the  residual  modes. 
Control  derivative  penalties  are  used  to  reduce  the  level  of 
excitation  of  these  modes.  Only  open  loop  maneuvers  between 
specified  states  in  a  specified  time  interval  were  considered 
and  the  results  indicate  that  the  problem  can  be  made  more 
practical  by  using  a  reduced  order  structural  model. 

In  [9]  Venkayyt  and  Tischler  report  on  the  effects  of 
modification  of  structural  parameters  on  the  dynamic  response 
of  flexible  structures  with  and  without  active  controls.  The 
control  problem  and  structural  optimisation  are  conducted 
independently  with  the  structural  optimization  based  on  the 
minimization  of  the  mass  with  the  fundamental  frequency  as  a 
constraint.  The  procedure  used  addresses  the  optimal  control 
design  with  a  quadratic  performance  index  involving  the 
energy  of  the  vibrating  structure  and  the  actuators. 
External  disturbances  which  initiate  or  continue  the 
vibration  of  the  system  are  not  considered.  A  state  space 
representation  for  the  open  loop  system  with  modal  reduction 
is  used  to  derive  an  expression  for  the  damping  ratio  and 
frequencies.  To  determine  the  optimum  feedback  control,  a 
quadratic  performance  index  with  specified  weighting  matrices 


is  minimized.  This  is  accomplished  through  the  solution  of 
the  steady  state  matrix  Riccati  equation. 


The  algorithm  used  for  the  structural  optimisation  is 
developed  in  detail  in  [6].  The  objective  is  to  minimize  the 
mass  of  a  vibrating  structure  subject  to  constraining  the 
fundamental  frequency  to  a  specified  value.  A  scaling  factor 
is  used  such  that  the  frequency  of  the  scaled  design  is  equal 
to  the  specified  fundamental  frequency  of  the  unsealed 
equation.  The  scaling  factor  is  iterated  from  0  to  1,  vith 
the  final  value  converging  to  the  desired  result.  The  paper 
shows  that  modification  of  the  structural  parameters  did  not 
significantly  alter  the  control  design  with  no  external 
disturbances,  but  the  changes  are  expected  by  the  authors  to 
have  a  more  pronounced  effect  when  an  external  disturbance  is 


Khot.  Venkayya.  and  Eastep  in  [3].  examine  the  dynamic 


stiffness  coefficients.  The  paper  uses  a  procedure  to  examine 
the  effect  of  structural  modifications  on  the  pointing  error 
of  a  tetrahedral  truss.  Essentially,  the  cross-sectional 


areas  of  the  tj 


optimized  to  mimimize  a 


displacement  function  and  in  a  separate  algorithm  are  used  to 
minimize  the  mass  of  the  structure,  while  constraining  the 
frequency  distribution.  An  iteration  process  is  used  until 
the  optimality  criterion  are  satisfied.  They  show  that  by 


changing  the  stiffness  of  the  members  based  on  the  optimum 


c 


design  criterion,  «n  inprovonent  in  the  dynaaic  response  of 
the  structure  to  a  external  disturbance  can  he  nade.  The 
iaprovenent  depends  on  the  criteria  need  and  the  initial 
conditions. 

In  [4],  Junkins,  Bodden,  and  lasat  develop  an 
eigenspace  optiaization  approach  for  design  of  feedback 
controllers  for  aaneuver s/ vibration  arrests  of  flexible 
structures.  They  state  that  aost  structural  and  control 
optiaality  criteria  can  be  stated  explicitly  or  iaplicitly  in 
teras  of  the  eigenvalues,  eigenvectors,  or  directly  as  a 
function  of  a  global  structural  and  oontrol  paraaeter  vector. 
The  eigenvalues  and  eigenvectors  are  shown  to  be  a  function 
of  the  structural  and  control  paraaeters  and  an  algoritha  is 
developed  to  coapute  their  partial  derivatives.  A 
differential  equation  aodel  for  the  Draper/RPL  Configuration 
is  used  to  deaonstrate  the  approach.  In  the  nuaerical 
exaaple  presented,  the  structural  paraaeters,  sensor 
locations,  and  actuator  locations  are  fixed  and  the  control 
gains  which  stabilize  the  first  few  aodes  with  the  saallest 
sua  square  are  solved  for.  A  performance  aeasure  and 
constraints  are  defined  in  terns  of  the  e i genso 1 ut i on  with 
explicit  dependence  upon  the  control  paraaeters.  The 
algoritha  is  applied  to  place  the  systea's  closed  loop 
eigenvalues  in  a  specified  region  snd  subject  to  this 
condition  ainiaize  a  robustness  aeasure,  the  sensitivity  of 
the  eigenvalues. 


( 


AX' 


Following  the  work  of  Joskins.  Bodden, snd  Kaset  in  [4] 
sn  sigsnspsce  optiaization  approach  is  used  to  incorporate 
optimal  control  into  the  strnotnral  design  process  of  large 
flexible  space  strnctnres.  An  expression  for  the  rates  of 
ekange  of  eigenwalnes  and  eigenvectors  with  respect  to 
parameters  of  the  system  is  required.  Tittriek  in  [9] 
determined  the  first  derivatives  of  the  eigenvalues  for  a 
self-adjoint  system  and  Fox  and  Kapoor  in  [10]  developed 
expressions  for  the  first  derivatives  ofthe  eigenvectors. 
Plant  and  Hnseyin  in  [7]  derived  general  expressions  for  the 
derivatives  of  eigenvalnes  and  eigenvectors  in  non-self- 
adjoint  systems.  Junkins.  et  al,  used  this  work  to  develop 
an  algorithm  to  minimise  an  optimality  criteria  and  satisfy 
constraints  stated  as  a  function  of  the  parameters  of  the 
system  (p)  and  the  eigenvalues.  As  stated  in  their  paper  a 
necessity  for  the  optimisation  algorithm  is  to  efficiently 
compute  the  partial  derivatives  of  the  eigenvalues  and 
eigenvectors  with  respect  to  the  structural  and  control 
design  parameters  (p). 

After  developing  their  algorithm,  they  applied  it  to  a 
mathematical  model  of  the  Draper/RPL  Configuration  and 
solved  for  the  optimum  gain  matrices  holding  all  the 
structural  parameters,  sensor  locations,  and  actuator 


locations  constant.  These  gains  sre  those  which  provide 
specified  values  for  the  first  few  model  frequencies  and 
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II.  PBmorifflSI  QE  MAlKMhUSAL  IfifiBL 

la  order  to  develop  an  analytical  aodel  it  is  accessary 
to  develop  sa  idealized  aodel  of  the  atraetaral  aystea  to  be 
atadied.  The  atractare  aaed  aa  a  aodel  for  this  report  ia 
the  'Draper  /  KPL  Coaf if aratioa*  aodel  illaatrated  ia  Figare 
1.  The  atractare  coaaiata  of  a  rigid  ceatral  hab  with  foar 
atractarally  ideatical  flexible  appeadagea  caatilevered 
radially  froa  the  hab.  The  ceatral  hab  ia  aaaaaed  to  be 
aapported  on  an  air  bearing  table.  Ezperiaeatal  work 
concerning  the  planar  rotational/  vibrational  dynaaica  of 
thia  U.8.A.F.  Socket  Propalaioa  Laboratory  dea on a t r a t i on 
aodel  ia  preaently  being  condneted  at  the  Charles  Stark 
Draper  Laboratory. [4] 

The  cantilevered  appendages  aodeled  here  can  be  thought 
of  as  a  continuous  beaa  with  a  luaped  aasa  at  the  end.  The 
differential  equations  of  notion  can  then  be  thought  of  aa 
the  natheaatical  language  which  describes  the  analytical 
aodel.  Thia  set  of  differential  equations  will  be  referred 
to  aa  the  aathenatieal  aodel  of  the  structure.  The 
natheaatical  aodel  developaent  begins  with  the  derivation  of 
the  equations  of  notion  for  an  unforced  aystea.  Expressions 
for  the  kinetic  and  potential  energy  are  developed  for  the 
aystea.  These  are  discretized  using  the  assuaed  nodes  aethod 
and  terns  of  order  higher  than  two  ignored.  The  aass  and 
stiffness  aatriees  are  then  identified.  Following  the 


I 

S' 


unforced  or  uncontrolled  aodel  i«  the  derivation  of  the 
controlled  systea.  The  control  gains  need  will  he  those 
fonnd  optiana  by  Junkins  in  [4].  The  next  step  in  the 
analysis  is  to  solve  the  differential  equations  to  obtain 
the  dynaaical  response  of  the  free  vibration  and  also  the 
response  of  the  controlled  systea. 

BBIlfllllgM  Of  IMJLS  FOR  EQUATIONS  QE  NOTION 

rj  is  a  position  vector  of  the  ith  appendage 
T  is  the  total  kinetic  energy  of  the  systea 
T^0  is  the  kinetic  energy  of  the  rigid  hub 
Tjjj  is  the  kinetic  energy  of  the  ith  appendage 
*  .  f.  i  represents  the  body  fixed  reference  fraae 
u.  v  are  the  deflections  of  the  appendages  in  the 
^  directions,  respectively. 

v  is  the  angular  rotation  of  the  body  fixed  reference 
fraae  with  respect  to  the  inertial  fraae. 

R  is  the  radius  of  the  hub 
L  is  the  length  of  the  appendages 
t  is  the  thickness  of  the  appendages 
h  is  the  height  of  the  appendages 

l20  is  the  rotational  aoaent  of  inertia  of  the  tip 

aa  as . 

aQ  is  the  aass  of  the  hub 

a*  is  the  aass  of  the  point  aasses 


» 


Aaummmtg 

Th •  loagitadinal  and  oat  of  plaaa  ribratioa  of  the 
appendages  ara  of  aaeh  higher  freqaeacy  than  tha  transverse 
▼ibratioaa  aad  ara  aaglaotad.  Aat i- sy aae tr i e  daforaatioaa 
vara  astaaad  aaeb  that  tha  daflaetioa  of  tha  firat  appaadaga 
vaa  aqaal  ia  aagaitnda  bat  oppoaita  ia  diraotioa  to  tha 
saeoad  appaadaga  aad  tha  daflaetioa  of  tha  third  appaadaga 
vai  aqaal  aad  oppoaita  to  tha  foarth,  aaa  figare  2.  Also 
since  tha  hab  is  sapportad  by  aa  air  baariag  tabla,  tha  forea 
of  gravity  was  aaglaotad  ia  tha  derivation  of  tha  potaatial 
aaargy.  Tha  daapiag  ooaffieiaat  vas  assaaad  to  ba  aqaal  to 
zero.  All  foar  appeadages  ara  aasaaad  to  ba  straetarally 
identical,  that  is 

L1  *  L2  “  L3  “  l4  “  L 

tl  -  t2  -  t3  -  t4  -  t 

kl  “  hj  ■  hj  ■  h4  ■  h 

aad  all  foar  tip  aassea  ara  idaatieal  and  aqaal  to  a2 . 


Figaro  2.  Aat i-Syaae trio  Deflections 


The  kinetic  energy  of  tke  hub  alone  is  just  that  of  a 


soliiLjBjl  inder  and  is  written  as 

Tb0  -  1/2  Ie  w*  (1) 

where  I0  «  1/2  »0 

To  deteraine  the  kinetic  energy  of  the  first  appendage 
a  position  vector  is  written  as 

-  x  ^  +  v  f 

The  tiae  rate  of  change  of  the  eleaents  of  appendage  1 
are  deterained  froa  the  first  derivative  with  respect  to  an 
inertial  reference  fraae.  All  derivations  here  are  with 
respect  to  an  inertial  reference  fraae. 

d  Fj/dt  ■  t  f  +  »  x  tj 

rj»-vw^  +  <  v  +  x  w  )  ^  (2) 

The  kinetic  energy  of  the  first  appendage  is  the  inner 

<1  JL 

product  or  dot  product  of  rj  .  r^  integrated  over  the  aass 
froa  the  end  of  the  hub  to  the  end  of  the  appendage. 


2  v  x  w  +  x2  w2  )  (3) 
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x-k+L 


Tbl  “  1/2 


r2  da 


x«R+L 


Let  Ij  *  I  xz  da 


Then  th«  kidetic  energy  for  tke  first  appendage  can  be 
written  as 


1/2  1,  w2  +  1/2 


(  v2w2  +  v2 


+  2  v  i  x  )  da 


The  kinetic  energy  of  tke  tip  aass  is  found  by  writing  tke 
position  vector  ,  taking  tke  appropriate  derivatives  and 
coaputing  tke  dot  product.  Tkerefore  : 


(  R  +  L  )  $  +  v(  x=R+L  )  ^ 


wkere  v  is  evaluated  at  x  *  R+L  Taking  tke  first 

derivative  yields 


v  V  + 


x  r2 


-  (  ▼  +  *  (R+L)  )  f  -  v  »  ^  (6) 

r2  *  r2  m  ^2 

*  v2  +  2  v  w  (R+L)  +  w2(R+L)2 

+  V2  w2 

The  angular  velocity  of  the  tip  •as*  is  the  angular  velocity 
of  the  body  fixed  reference  frame  plus  a  component  due  to  the 
deflection  of  the  end  of  the  appendage  where  the  tip  mass  is 
located. 


where 


r2  -  w  +  v  » 


v  '  »  d  (  8v  /  dx  )  /  dt  evaluated  at  x  ■  R+L  (7) 


The  kinetic  energy  of  the  tip  mass  can  now  be  expressed  as 


T„j  -  1/2  .j  Tj2  +  Il0  ,22 


While  the  rotational  moment  of  inertia  of  a  tip  mass 
is  small ,  it  will  be  included  in  the  general  formulation. 

In  a  similar  manner,  the  kinetic  energy  of  the  third 
appendage  and  its  tip  mass  were  derived  and  are  shown  here  as 


y-R+L 


Tb3  -  1/2  I3  w2  +  1/2  /  ( u2 w2  +  u2 


+  2u  w  y)  dm  (9) 


where 


y=  R+L 


and  the  expression  for  the  tip  mass  of  the  third  appendage  is 

TH  -  1/2  H  V  ♦  J2c  V 

where 

V42  =  u2  +  2  u  w  (R+L)  +  u2  w2  +  w2(R+L)2  (10) 


and  w4  *  w  +  u  * 

where  i  '  *  d  (  3  u  /  Jy  )  /  dt  evaluated  at  y  -  R+L 


Now,  since  only  the  case  of  antisymmetric  deformation 
is  being  considered  ,  where 

n3  =  -u4  ,  y1  -  -v2 

it  follows  that  the  kinetic  energy  of  the  first  appendage  and 
tip  mass  is  eqnal  in  magnitude  to  the  kinetic  energy  of  the 
second  appendage  and  its  tip  mass.  Likewise  the  kinetic 
energy  of  the  third  appendage  and  tip  mass  is  equal  in 
magnitude  to  the  kinetic  energy  of  the  fourth  appendage  and 
point  mass.  Therefore,  the  total  kinetic  energy  of  the 


system  can  be  written  as  the  sum  of  the  individual  terms  as 


T  »  Tb0  +  2  Tbl  +  2  Tb2  +  2  Tb3  +  2  Tb4 


To  expand  this  expression  let 


Ic  ■  1/2  I#  +  2  Ij  +  2  m2  (R+L)  2 


The  total  kinetic  energy  can  be  regrouped  as 


T  *  (Ic+2  I2c)w2  +  pt  h  /  (u2w2+u2+2n  w  x)  dy 


/ 


/ 


+  p  t  h  /  (v2w2+v2+2  v  w  y  ) dx 


+  m2 (  u2  +  2  u  w  (R+L)  +  u2w2)  +  I2c(2wu*  +  u'2 


(11) 


+  m2(v2  +  2  v  w(R+L)  +  v2w2)  +  I2c(2wv'  +  v'2) 


Similar  to  the  work  of  Turner  and  Chun  in  [11]  a 


0  <  s  <  L 


Likewise  for  eppeadeges  3  and  4  , 


v(t.z)  =  £  Vi(t)  *(z) 


(15) 


The  assumed  modes  method  leads  to  a  generalized  parameter 
model,  where  the  U^Ct)  and  Vj(t)  are  the  generalized 
displacement  coordinates.  Taking  the  first  derivative  of 
equations  (14)  and  (IS)  with  respect  to  time  yields 


P  . 

u(t,z)  -  V  U^t)  +i(z) 
1-1 


v(t,z)  =■  £  Vj(t)  ^i(z) 

i=l 


Squaring  the  above  equations 


P  p 

u2(t,z)  «  2  ^  ®i  +i 


v2 ( t ,  z)  .  £  £  v,  y3  ^ 

i«i  j-i  J  J 


The  above  expressions  for  the  generalized  coordii 


(16) 


(18) 


(19) 


their  tii 


tes  of  change,  equations  (14), (15),  (16).  (17), 


(18),  and  (19)  can  be  inserted  into  the  expression  for  the 


it  t  f . 


total  kinetic  energy  of  the  systea,  equation  (12). 


i 

i 


I 


Potential  Energy 

The  potential  energy  of  the  ayaten  can  be  written  aa 
the  ana  of  the  grayitional  and  the  atrain  energy  or  elaatic 
potential  energy. 


V  *  *  Vtl.,tlc  (22) 

• 

Since  the  hnb  ia  aaanaed  to  be  anpported  by  an  air  bearing 
table  (  or  in  orbit  )  the  gravitional  potential  will  be 
C  neglected  and  the  potential  energy  written  in  teraa  of  the 

atrain  energy  alone. 


c 


fc- 


v  ■  <”> 

The  total  atrain  energy  ia  the  ana  of  the  contribntiona  froa 
all  four  appendages.  The  expression  for  veiagtic  has  been 
derived  extensively  in  the  literature  aa 

L 

V.L.ti.  *  2/2  E  I  j  (  >* 

0 

or  for  all  four  appendages  (  1*2,  3*4  ): 


I 


velastic  E  1 


# 


C 


(22) 


where  u"  end  ▼"  ere  the  second  partial  derivatives  of  n  and 
v  with  respect  to  z,  respectively.  E  is  the  modulus  of 
elasticity  of  the  arms  ,  also  known  as  Young's  modulus,  and 
is  a  property  of  the  material.  I  it  the  area  moment  of 
inertia,  based  only  upon  t  and  h  as 

I  -  1/12  h  t3 

Again,  the  assumed  modes  method  leads  to  the  generalized 
parameter  model  with  U(t)  and  V(t)  aa  the  generalized 
displacement  coordinates  and  ^^(z)  as  the  admissable 
functions. 

u' '  and  v  *  *  he  come 

P 

u"  «  £  Dj  ( t )  (z)  (23) 

i-1 


P 

V ’  -  £  Vi(t>  ^i' '  (*)  (24) 

i-1 

where  ia  the  second  partial  derivative  of  ^  with 


s* 

s‘ 

i 


respect  to  z,  see  equation  (13). 

♦  t"  -  (  iji/L) 2  cos  (  iits/L)  +  (-1)  i+1(in/L)2  (25) 

Utilising  equations  (23)  and  (24)  and  after  grouping  terns 
and  arranging  in  matrix  notation  ,  equation  (22)  becomes 

^elastic 


(26) 

Recalling  Lagrange's  equation  for  the  derivation  of  the 
differential  equations  of  motion  as 

d  OL  /  d  q  j  )  /  dt  -  d  L  /  dqj  -  Q 

i"l ,  .  .  .  no.  of  generalised  coordinates 

where  L  ■  T  -  y  and  Q  represents  the  noncon s e rv a t iv e 


forces. 


Applying  this  aethod  to  the  expressions  for  the  kinetic 
energy  end  potential  energy  of  the  free  vibration  system  end 
negleeting  the  terms  of  order  three  or  higher  and  writing  the 
kinetic  and  potential  energies  in  the  fora  : 

T  -  1/2  ?T  N  t 

V  »  1/2  qT  K  q 

where  M  and  K  are  called  the  mass  and  stiffness  matrices, 

T 

and  q  *  [  0,  U2  .  .  .  Up,  V2  .  .  .  Vp  ] 

The  elements  of  the  mass  matrix  are  given  by  : 


N  (1.1)  -  2  (Ic  +  I2c  ) 


E+L 

/ 


M ( 1 ,  i+1 )  ■  2  pt  h  /  x  t  dx  +2  (E+L)  a2 

E 

+  2  *2c  ♦,<*-*+L) 


XI,  i+p+1) 


2  p  t  h 


E+L 

f  y  ♦  dy 

E 


+  2  m2(E+L)  ♦(y«E+L) 


+  2  I2c  ♦»(y-E+L) 

M(  i+1,1)  -  N( 1 , i+1) 

H( i+p+1 , 1 )  -  1(1, i+p+1) 


E+L 

(♦. 


dx 


M(i+l,j+l)  -  2  p  t  h 


+2  a2  ^  j  (  x « E+  L )  ♦  jtR+L) 


♦ 2  I2c  (♦i*(x-R+L))  (^Mx-R+D) 


M(  i+p+1 , j+p+1 )  «  2  p  t  h 


X+L 

J  *i  *)  *7 


m2  ♦  jU-R+L)  ♦j(x-R+L)  +  2  I2c  +  j  '  ( x-R+Ll  f  j  » ( x«R+L) 


M( i+1. j+p+1)  -  M( i+p+1 . j+1)  «  0 


(27) 


where  i*l,...  p  end  j-1,...  p 


The  resnlting  axes  aatrix  is  syaaetrie  end  of  order  2  p+  1 
where  p  is  the  nnaber  of  aodes  being  ewalneted.  The  eleaenti 
of  the  stiffness  aetrix  follow  froa  equation  (2d)  as 


K(l,l)  -  1(1, j+1)  -  I(l.j+l+p)  -  0 
K( i+1 , 1)  -  K( i+l+p , 1)  -  0 


1( i+1 , j+1 )  -  2  E  I 


L 

/  ♦i’<* 


)  ♦j'Mx)  da 


I( i+l+p. j+l+p)  -  2  E  I  I  ♦£”(*)  ** 


•MvV. 


i-:,. t> 


:  *>*  - 


I<i+l.j*l+p)  -  Ki+l+p.j+l)  -  0 


(28) 


Tl*  victor  differential  aquation  governing  the  motion 
of  the  naoontrolled  tpctca  takes  the  familiar  form 

M  q  +  K  q  -  0  (29) 

where 

4  -  l  9  .  Dj  .  .  .  Up.  Vx  .  .  .  Vp  JT 

note:  by  assnming  anti-aymme trie  motion  of  the  appendages  the 
number  of  generalised  coordinates  is  redneed  by  2  p  . 

The  nest  step  in  the  analysis  is  to  solve  the 
differential  equation  (29)  to  obtain  the  response  of  the 
uncontrolled  system.  The  method  need  was  to  solve  the 
familiar  eigenvalue  problem, 

K  x  -  X  N  i  (30) 

for  the  eigenvalues  and  eigenvectors  for  the  system.  The 
first  10  modes  were  analysed  and  the  resulting  eigenvalues 
and  eigenvectors  are  listed  in  Appendix  A  along  with  a  plot 
of  the  first  seven  of  those  modes  depicting  the  relative 
motion  of  the  free  arms.  The  modes  appear  ss  pairs,  very 
close  in  frequency,  depicting  what  Junkins  in  [4]  calls 
'opposition'  and  'unison'  modes.  The  opposition  modes  are 


siaple  cist il «v*r  btta  aodti  ehariottrit «d  by  tk«  idj  teaat 
betas  loving  in  opposition.  Ho  ititis  the  nnison  nodes  its 
perturbed  enntilevor  nodes,  with  jnst  slightly  higher 
frequencies,  with  ell  four  betas  nosing  in  unison  tnd  the  hub 
hewing  non-tero  rotition.  The  rotition  of  the  hub  it  to 
counter  the  unison  aoveaent  of  the  tppendtges  tnd  point 
atsses  to  conserve  tngultr  aoaentua  in  the  systea.  The 
nuaerictl  vtlues  of  the  frequencies  for  the  first  seven  nodes 
ere  very  close  to  those  listed  by  Juakins  in  [4]. 


Figure  4.  Opposition  tnd  Unison  Nodes 


gflMIIfllLM  5LSXBM 

The  equation  oi  otion  for  the  controlled  linear  systea 
can  be  written  in  tbe  general  fora  of 

II  q  +  C  t  +  I  4  -  B  1  (31) 

where 

N  ie  the  N  x  N  (N»2p  +  1)  eyaaetric  and  positive 

definite  aass  aatrix  derived  in  the  previous  section. 
C  is  the  N  x  N  structural  daaping  aatrix. 

K  is  the  N  x  N  syaaetric  positive  s ea  i-de f i ni t e 
stiffness  aatrix. 

B  is  an  N  x  a  control  aatrix. 

q  is  the  vector,  diaens i onN,  of  generalized 
coordinates  derived  froa  the  assuaed  aodes  aethod 
of  the  previous  section. 

M  is  the  control  vector  of  diaension  a. 

Only  the  case  of  linear  output  feedback  control,  with  local 
position  and  velocity  aeasureaents  available,  will  be 
considered  here.  In  this  case  the  gain  aatrices  are  constant 
such  that 

B  ■  -Oj  Hj  q  -  02  H2  q  (32) 

where  I|  and  H2  are  the  sensor  aatrioes  that  representthe 
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linear  relationship  of  the  locally  Measured  position  and 
velocity.  Snbstitnting  equation  (32)  into  equation  (31) 
yields  the  closed  loop  systea  equation 


i  *  g 

M  q  +  C  q  +  K  q 


(33) 


where 


C*  -  C  +  B  G2  H2 


or  since  for  the  problea  ezaained  zero  syatea  daaping  has 
been  assuaed: 


and 


C*  «  B  G2  H2 


I*  *  I  +  B  Gj  Hj 


The  next  step  is  to  solve  the  differential  equation 
(33)  governing  the  controlled  aystea  to  obtain  the  dynaaical 
response.  Equation  (33)  is  a  second  order  aatrix 
differential  equation  representing  the  closed  loop  systea. 
This  second  order  differential  equation  nay  be  written  in 
first  order  fora  as 


A  » 


D  y 


where 


-M 


(34) 


I*  C* 
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This  first  order  differeatisl  equation  can  be  solved  by 


assuaiag  a  solation  of  the  fora 


f 


(35) 


and  solving  for  the  eigenvalues  and  right  and  left 
eigenvectors  of  the  controlled  systen.  Since  aatriz  A  and  D 
are  of  order  2  N  ,  there  will  also  be  2  N  eigenvalues  and 
the  right  and  left  eigenvectors  will  be  of  order  2N.  They 
are  found  froa 


right  A  0|  *  D  0 | 


left  :  Aj  AT  ■  DT  Ti 


where  i-1,  .  .  .  2N 


(36) 


The  solution  to  the  uncontrolled  system  was  found 
using  the  values  for  6^  and  G2  found  optiaua  by  Junkins  in 
[4]  and  using  the  saae  expressions  for  B,  Hj,  and  Hj • 


8.60 

-1 .08 

-3  .09 

-1  .08 

-3  .09 

- 

21.80 

-1  .74 

-4.94 

-1  .74 

-4.94 

1  21 .90 

-1.74 

-4.94 

-1.74 

-4.94 

m 

-.238 

-  .016 

-  .050 

-  .016 

-  .050 

« 

.643 

.233 

.119 

-  .015 

.028 

.643 

-.015 

.028 

.233 

.119 

Following  the  seme  configuration  of  the  Draper  model  as 
Junkins.  torqne  actuators  were  admitted  on  the  central  hnb 
and  at  some  position  (z)  on  each  appendage.  The  actuators 


consisted  of  a  torqne  nj  applied  to  the  hnb.  a  torqne  at 
z  =  L/2  for  arms  1  and  2  and  a  torqne  jij  applied  at  z  =  L/2  for 


third  and 

four  th 

appendages . 

Matrix 

B  is 

then 

1 

2 

2 

j 

**1 

B  - 

0 

2  ♦*(z«L/2) 

0 

P2 

0 

0 

2  *'(z« 

L/2) 

•*3 

(37) 


where  +'(z)  as  stated  before  is  d(+(z))/dz  and  the  number  of 
elements  of  i«  equal  to  the  number  of  modes  being 
examined.  B  is  an  N  x  3  matrix. 

Co-located  position  and  velocity  sensors  are  located  at 
the  same  positions  for  appendages  1  and  2  and  appendages  3 
and  4.  This  has  the  effect  of  rendering  “  H2  *  H* 


1 

H 

lo 

0T 

0 

♦  T(z-L/2) 

0T 

0 

♦  T(z-L) 

oT 

0 

oi 

H 

♦T(z=L/2) 

0 

Ol 

H 

*T(z=L) 

(37) 


hence  H  is  a  5  x  N  matrix. 

Inserting  the  expressions  for  6^  , 
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,  v. 


*•  I  fairtniV  ■ 


G2  ,  B  and  H  into 


Hub  radius,  R 

1  ft  9 

Rotary  Inertia  of  Hub 

8  slug-ftz 

Mass  density  of  beams, p 

5.22  slug-ft3 

Modulus  of  Elasticity, 

E 

1.584E+09  lb/ft2 

Appendage  thickness,  t 

.0104166  ft 

Appendage  height,  h 

.5  ft 

Appendage  length,  L 

4.0  ft 

Tip  mass 

.156941  slug 

Rotary  Inertia  of  tip  mass 

.0018  slug-ft2 

Table  1.  Configuration 

Par  ame  ter  s 

for  Draper  Model 

The  parameter  optimization  technique  or  method  used  to 
minimize  the  mass  of  the  'Draper  /RPL  Configuration*  model 
involves  the  eigenspace  optimization  approach  developed  by 
Junkins  in  [4]  and  is  based  upon  the  work  by  Plant  and 
Hnseyin  in  [7].  As  in  any  optimization  scheme  ,  an  objective 
function  is  sought  to  be  minimized  subject  to  certain  imposed 
equality  and  inequality  constraints.  'An  algorithm  that  can 
effectively  optimize  a  global  structural  and  control 
parameter  vector  (  over  some  admissible  set,  to  minimize  some 
optimality  criteria  and  satisfy  constraints  stated  as 
functions  of  that  vector  and  the  eigenvalues  )  provides  a 
direct  method  for  c o n t r o 1 1 e r / s t r u c t u r e  optimization 
problems. * [4] 

SOLUTION  PROCEDURE 

Let  d  be  a  vector  of  global  structural  and  control 
parameters  such  that 

M  *  M(d)  ,  C*  =  C* ( d )  ,  K*  =  K*(d)  (39) 


it  follows  that  the  eigenvalues  and  eigenvectors  will  also  be 
a  function  of  d  . 
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The  eigenvalues  of  *  sy stea  can  be  a  highly  nonlinear 
function  of  the  structural  and  control  parameters  and  may 


contain  bifurcation  points  at  repeated  roots.  Theref or e, the 
constraints  should  include  a  measure  of  eigenvalue 
sensitivity  to  prevent  encountering  these  singular  points. 
Another  problea  described  by  Junkins  in  [4]  that  can  be 
encountered  is  in  attempting  to  place  eigenvalues  at  exact 
postions  or  'pole  placement'.  Therefore,  instead  of 
attempting  to  prescribe  an  exact  location  for  an  eigenvalue, 
a  region  (in  the  left  half  plane)  will  be  specified. 

The  technique  used  in  the  optimization  is  based  upon 
the  necessity  to  effectively  compute  the  eigenvalue 
sensitivities  to  changes  in  the  structural  parameters,  or 
better  stated  as  the  partial  derivatives  of  the  eigenvalues 
denoted  below  by: 

d\i  /  adx  i-1  .  .  .  N  (41) 

1=1  .  .  .  no.  of  parameters 

for  the  generalized  eigenvalue  problem  of  equation  (36). 

Plant  and  Huseyin  in  [7]  and  Rogers  in  [8]  present  an 
excellent  derivation  of  equation  (41),  the  results  of  which 
are  presented  below 

aXi/ad!  =  ?iT  [dD/adj  -  kjaA/adj  ]  (42) 


The  sensitivities  of  the  eigenvalues  can  be  computed  once  the 


eigenvalues  and  eigenvectors  have  been  found.  Although  the 
calculations  do  not  appear  to  be  eztreaely  difficult,  they 
are  tedious  and  tedious  work  is  best  left  to  a  aachine,  in 
this  ease  the  Cyber  computer. 

Expanding  equation  (42)  yields  : 

dkj/3dj  my  3D/ddj  a  ^  A/3  djo  j  (43) 

Remembering  that 


D  - 


Then 


3D/3d1 


3M/3d1 


[aK'/ddj  dcVadj 


3  A/3  dj 


-dM/ddi  0 


0  -3N/3d1 


(44) 


expanding  the  above  terms  and  using  the  chain  rule 
tensively,  the  partial  derivatives  needed  in  (44)  can  be 


calculated. 


K  -  1(d)  +  B(d)  Gxtd)  H(d) 


dK*/dd1  -  31/ddj  +  dB/ddj  G1  H  +  B  aGl/ddj  H 


+  B  G1  dH/ddj 


(45) 


i&d  C  *  B  62  H 


dC*/9d1  -  dB/dd|  G2  H  +  B  dG2/3d1  H  +  B  G2  dE/dd1 


(46) 


Once  the  initial  strnctnral  and  oontrol  pi 


specified,  the  sail  and  stiffness  aatrices  can  be  coapnted 


and  fed  to  an  algorithm  for  computing  the  eigenvalues  and 


eigenvectors  (left  and  right).  With  a  few  additional 


moderate  calculations,  the  sensitivities  of  the  eigenvalues 


can  also  be  computed.  Therefore,  the  requirement  to 


effectively  calculate  the  eigenvalue  sensitivities  has  been 


fulfilled  and  the  task  of  accomplishing  the  stated  objective 


can  be  tackled. 


The  objective  is  to  minimise  the  mass  of  the  structure 


given  an  initial  design  with  an  optimum  control  system 


already  developed  for  the  initial  design,  by  optimising  the 


structural  parameters.  The  objective  function  is  subject  to 


the  constraints  of  maintaining  specified  natural  frequencies 


and  damping  factors  eonstant  (or  within  specified 


tolerances).  The  mass  of  the  hub  will  be  held  constant  as 


will  the  tip  mass  at  the  end  of  each  appendage,  such  that 


k 

|r 


B 


only  the  aass  of  the  appendages  will  be  ainiaiaed.  Restated, 
the  objective  ie  to  ainiaize 

4 

M  -  P  ti  hj  Lj  (47) 

i»l 

subject  to  specified  natural  frequencies  and  daaping  factors. 
The  natural  frequencies  and  daaping  factors  are  related  to 
the  structural  paraaeters  through  the  eigenvalues  as 

«ni  ■  [  (Real  X*  )2  +  (Iaag  XA)2  J1/2 

*  -Real  Xj  /  «ni 
where  Xj  -  Xt(d) 

(48) 

The  role  of  the  eigenvalue  sensitivities  becoaes 
readily  apparent  when  the  gradients  or  partial  derivatives  of 
the  natural  frequencies  and  daaping  factors  with  respect  to 
the  structural  paraaeters  are  calculated  using  the  chain 
rul  e . 

dw^/ddj  ■  [dReal  X^/ddj  +  J IaigX j/J dj ] /  wni 

35/adj  -  - [d  Real  X^d^  /  *ni 

+  [  dReal  Xi/dd1  +  d Iaag  Xi/ddl  J  /  **ni3/2 

(49) 

If  the  structural  paraaeters  are  considered  to  be  the 


thickness,  t,  the  height,  h,  end  the  length,  L,  of  each 
appendage  then  the  paraaeter  vector  d  becoses 


^  B  t  f  j  •••  tg  ,  hj  •  •  •  hg ,  Lj  •  •  •  Lg  ] 


(SO) 


In  the  problea  considered  here,  all  four  appendages  are 
kept  structurally  identical,  such  that  the  paraaeter  vector 
is  reduced  to  only  three  eleaents 


d  -  [  t  h  L  ] 


(SI) 


If  a  need  were  to  arise  to  change  certain  appendages 
while  holding  others  constant  or  to  ohange  all  four,  but 
weighted  differently,  then  the  diaension  of  d  would  increase 
accordingly.  Referring  back  to  equations  (45)  and  (46)  the 
gain  aatrices  Gj  and  62»  functions  of  fhe  control  gains 
alone  and  therefore  the  partial  derivatives  of  G2  and  G2  with 
respect  to  the  current  paraaeter  vector  are  sero.  That  is 


flGj/gd  «  0 


«G2/8d  -  0 


(52) 


Also  the  control  aatrices  B  and  H  are  not  functions  of  t  or 
h,  therefore 


dB/dt  -  dB/dh  -  0 


dE/dt  -  dH/3h  -  0 

(53) 

This  reduces  the  nsiber  of  calculations  required  when 
computing  the  eigenvalue  sensitivities  and  hence  the 
gradients  of  the  natural  frequencies  and  damping  factors. 

Therefore,  the  partial  derivatives  of  aatrices  A  and  D 
reduce  to  calculating  the  partial  derivatives  of  the  mass  and 
stiffness  aatrices  for  eleaents  t  and  h  of  the  paraaeter 
vector  d. 

For  eleaent  1.  the  thickness  t  : 


9M(l.l)/at  -  4/3  p  h  Ij  (I1  is  defined  in  ch.2) 

ana, j+D/at  -  2  P  h  pXj 

3M(i+j,i)/at  -  ' ' 

an( i+j+p, 1 ) /at  - 


vhere  j*l,p  and  p  is  the  nuaber  of  aodes. 


and 


P*j 


R+L 


da 


R 

dM(i+l,j  +  D/  dt  -  2  p  h  Pj  j 

a«( i+i+p, j+i+p)/at  « 


vhere  i*l,p  and  j"l,p 


(54) 
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9 

L 

f 

9 

and  pi.j  -  I  *i<z)  *j(*>  d* 

0 

a 

9 

Similarly  for  the  stiffness 

r 

t  r  O  .  r, 

s“ 

M 

i 

*  0 

3K(  i+1 ,  j  +  1 ) /3  t  =  .5  E  h  t2  Pifj" 

u 

■i 

e 

* 

9 

3K(i+l+p. j+l+p)/3t  « 

i 

Si 

1 

L 

r 

I 

\ 

where  pi,j"=  /  *j"<z)  dz 

t  1 

c 

. 

0 

j 

* 

(56)  < 

t 

9 

For  eleaent  2,  the  height,  h  : 

■ 

* 

4 

9 

dM(l,l)/ah  -  4/3  P  t  Ij 

3M(l.j+l)/3h  -  2  p  t  P*j 

3M( 1 , j+l+p) /3  h  - 

dM(l+j,l)/dh  - 

** 

. 

• 

* 

> 

* 

« 

«* 

3M< j+l+p, l)/3h  -  • » 

i 

J-l.P 

( 57 ) 

3M( i+1, j+l)/3h  -  2  p  t  Pk  j 

► 

*> 

► 

!C 

■ 

3M( i+l+p,j+l+p)/3h  -  * • 

£ 

(58)  | 

. 

ell  other  elements  of  M  ere  zero. 

» 

a* 

•  j 

*' 

For  the  stiffness 

fc 
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y.\v ,% -,v/.  /. ,%  /•  .*•»  v>  \  -  •>  V  v  v  v  v  *,*  r„* 

c 

.*  :■  ^  v-y-vfvCv  n-  wo  .i-.-'s, 

8I( i+1. j+l)/8h  -  1/6  E  t3  P 


8K( i+l+p, j+l+p)/8h  -  *  * 

i-l.p  .  j-l»P 

and  all  other  elements  of  K  are  zero. 


(59) 


For  the  third  eleaent,  L  ,  the  calcnlationa  become 
considerably  sore  involved  as  B  and  H  and  4  are  functions 
of  z  which  ranges  from  0  to  L. 


dD/dL  - 


8A/3L  = 


f: 


8M/8L 

8C*/8L 


-8M/8L 

0 


-8M/8L 


(60) 


Expanding. 

8N(1,1)/8L  -  4p  t  h  (R+L)2  +  8  n2(R+L) 
8M(l,j+l)/8L  *  2  p  t  h  8PXj/3L 

+  2  n2  4j(L)  +  2  I2c  8  ♦jMD/JL 
8M ( 1 , j  +  l+p ) / 8L  =  *  ' 

8M( j  +  l.l)/3L  «  '  ' 

8M(j+p+l,l)/8L  -  " 

j“l»P 


and  8fj'(L)/dL  «  -<-l)j+1  (jn/L)2 

8M( i+1, j+l)/3L  -  2p  t  h  8P .  ,/8L 

1  »  J 

+  2  I2<J  8(4i'(L)  f  j  '  (  L)  )  /  3  L 


dM<i+l+p, j+l+p)/3L 


9 


where  3  <  ^  '  +j’)/3L  =  -2  (in  jn)2  (-1) i+1  (-1)J  +  1  /L3 

(61) 


Expending  the  remaining  terse  in  aetrix  D  yields  : 


3I*/3L  -  3I/3L  +  3B/3L  Gj  H  +  B  Gj  3H/3L 
9C*/3L  «  3B/3L  G2  B  +  B  Gj  3H/9L 

3K(  i+l,j  +  l)/3L  -  1/6  E  h  t3 
3I( i+l+p. j+l+p)/3L  -  " 

where 


P 


i.  j 


9  9  m 


L 

J  dz 

0 


3B(  i+1 ,  2 ) /3L  -  2  3$*  •  (  z-L/  2) /3L 
9B(  1+1+p,  3  )  /3L  - 

where 

dfj* ( z-L/2)/3L  -  -(in/L2)  sin(in/2) 


3 H(  1 ,  j  +  1 )  / 3 L  -  3#j  (  z-L/2)/3L  -  0 
3H(4. j+l+p)/3L  - 
3  H(  3  ,  j  +  1 )  /3  L  -  3^j(z-L)/3L  -  0 
3H( 5  » j+l+p) /3L  - 


(62) 


(63) 


The  above  equations,  although  lengthy  when  expanded,  involve 
repitition  and  matrix  multiplication  which  are  ideally  suited 


for  a  computer. 


PI  COMPUTER  PROGRAMS 


Utilizing  the  above  developaent,  a  computer  program  was 
written  to  calculate  the  gradienta  of  the  natural  frequency 
and  the  damping  factor  for  each  eigenvalue.  The  program  is 
listed  in  Appendix  C  .  It  begins  by  calculating  the  mass 
and  stiffness  matrices  for  a  given  d  ,  initial  values  of  t. 

h,  and  L.  The  control  matrices.  B  and  H,  are  computed  and 

« 

combined  with  the  given  gain  matrices  to  form  the  A  and  D 
matrices  of 

D  a  =  X  A  a 

and 

DT  y  -  A  AT  r  (50) 

to  compute  the  eigenvalues  and  left  and  right  eigenvectors. 
Once  A  and  D  are  formed  they  are  input  to  an  IMSL  subroutine 
called  'EIGZF'  along  with  the  appropriate  dimensions.  EIGZF 
returns  the  eigenvalues  and  corresponding  eigenvectors.  The 
transpose  of  A  and  D  are  computed  and  also  input  to  EIGZF, 
which  returns  the  left  eigenvectors. 

Subroutine  CALCBA  is  called  to  compute  the  gradients  of 
the  A  and  D  matrices  with  respect  to  the  input  structural 
parameters.  This  information  along  with  the  eigenvalues  and 
eigenvectors  is  passed  to  another  subroutine  called  GRADNT, 
which  as  the  name  implies,  computes  the  gradients  or 
sensitivities  of  the  natural  frequency  and  damping  factors. 
The  eigenvalue  sensitivities  are  calculated  by  expanding 


equation  (42)  and  using  the  input  values  for  the  eigenvalues, 
left  and  right  eigenvector  a,  and  the  gradients  of  aatrices  A 
and  D.  The  partial  derivatives  of  the  natural  frequencies 
and  daaping  factors  can  nov  be  caluculated  using  equation 
(49)  . 

There  are  aany  standard  software  optiaization 
algorithas  available  that  allow  the  user  to  ainiaize  a  stated 
objective  function  subject  to  specified  equality  and 
inequality  constraints.  The  algoritha  used  here  was  one 
called  CONNIN  (constrained  ainia iz ation)  which  is  described 
in  [12].  It  is  a  FORTRAN  prograa  in  subroutine  fora  for 
constrained  function  a inia iz ation.  The  user  aust  provide  a 
aain  prograa  to  call  CONMIN  and  evaluate  the  objective  and 
constraint  functions  along  with  the  associated  gradients. 
CONMIN  then  uses  this  inforaation  in  conjunction  with  the 
aethod  of  Feasible  Directions  to  solve  the  constrained 
probl ea. 

A  aain  prograa  was  written  to  call  CONMIN  and  perfora 
the  necessary  calculations.  The  prograa  initializes  all 
variables  coanon  to  CONMIN  and  then  proapts  the  user  for 
certain  inputs.  Those  inputs  are  DALFDN,  DABFUN,  CT,  THETA 
which  specify  the  accuracy,  tolerances  ,  and  how  far  off  to 
aove  froa  the  initial  values  of  the  decision  variables  (in 
this  case  the  decision  variables  are  the  structural 
paraaeters).  Next  the  user  is  asked  to  input  the  initial 
values  for  the  decision  variables.  Subroutine  CALOBJ  is 


called  which  calculates  the  objective  function 

OBJ  =  4  p  t  h  L  (64) 

where  p  is  the  mass  per  unit  volume 

Subroutine  CALOBJ  calls  the  program  described  earlier  to 
provide  the  eigenvalues  and  gradients  of  the  natural 
frequencies  and  damping  factors  (which  are  necessary  since 
the  equality  constraints  are  the  natural  frequency  and 
damping  factor  and  the  objective  is  to  minimize  the  mass 
while  holding  them  constant).  The  constraints  then  become 

G  ( 1 )  =  oi  j  -  o>Q 

® ( 2 )  *  5i  -  gD  (65) 

where  o>j  and  gj  are  the  original  values  of  the  natural 
frequency  and  damping  factor  calculated  for  the  initial 
structural  parameters  and  the  gains  found  optimum  by  Junhins 
in  [4].  These  are  the  values  desired  to  be  held  constant 
while  £0  and  u>Q  are  the  values  re-computed  during  each 
iteration  of  CONNIN  and  evaluated.  The  gradients  of  the 
natural  frequency  and  damping  factor  are  also  passed  to 
CONNIN  to  provide  convergence  information.  The  constraint 
gradients  provide  a  direction  for  the  process  such  that  the 
value  of  the  objective  function  continues  to  decrease  while 
satisfying  the  constraints  G(l)  and  G(2).  The  ideal 
performance  of  obtaining  G(l)  or  G(2)  equal  to  zero  can 
seldom  be  obtained  in  a  digital  computer,  therefore  CONNIN 


uses  a  constraint  thickness  parameter  to  define  the  region 
over  vhich  G(l)  or  G(2)  can  be  assumed  zero  for  optimization 
purposes.  This  is  one  of  the  prompted  input  values  called 
CT,  so  the  program  allows  the  user  to  input  how  stringent  the 
process  is  to  be.  The  constraints,  G,  are  considered  active 
if  they  are 

CT.  LE  .  G( j ) .  LE  .  ABS  ( CT ) 

where  CT  is  the  constraint  thickness  parameter.  If  they  are 
no  longer  active,  the  specified  upper  and  lower  bounds  on  the 
decision  variables  (structural  parameters)  may  become  active. 
When  the  objective  function  can  no  longer  decrease,  the 
process  is  terminated. 

Therefore,  the  result  should  be  a  structure  that  begins 
with  a  specified  feasible  set  of  structural  parameters, 
optimizes  the  control  system  ,  and  then  optimizes  the 
structural  design.  This  procedure  can  be  repeated  until  the 
desired  mission  objectives  are  achieved. 


IV.  ANALYSIS  Oj?  RESULTS 

The  objective  of  this  report  is  to  minimize  the  mass  of 
the  appendages  while  constraining  either  the  fundamental 
frequency  or  another  natural  frequency  to  predetermined 
values.  The  gradients  of  the  eigenvalues  were  calculated, 
with  respect  to  each  element  of  d  ,  for  the  firstten 
eigenvalues.  This  information  was  used  while  examining  four 
different  cases.  Those  cases  are:  constraining  the 

fundamental  frequency  to  the  initial  value  found  using 
J unkins'  optimum  control  law,  constrain  the  third  natural 
frequency  to  an  initial  value,  constrain  the  third  natural 
frequency  to  an  acceptable  region,  and  constraining  the  third 
natural  frequency  to  a  specified  region  while  also 
constraining  the  fundamental  frequency. 

GRADIENTS  AND  EXPECTATIONS 

The  program  used  to  call  CONNIN  was  modified  to  print 
the  gradients  of  the  eigenvalues  and  was  run  separately.  The 
gradient  of  the  first  eigenvalue  for  all  three  parameters  was 
found  to  be  approximately  -2.0  E-10.  This  would  imply  that 
the  'surface'  is  very  flat  and  convergence  to  a  new  minima 
difficult.  Therefore,  one  would  expect  very  little  change  in 
the  objective  function  for  an  input  of  initial  values  and 
difficulty  in  convergence  if  a  non-feasible  solution  was 
input . 

The  second  eigenvalue  is  the  conjugate  of  the  first,  so 


the  gradient  of  the  third  eigenTtlte  was  eztaiaed.  The 
gradient  with  respect  to  the  thickness  was  fonnd  to  be  very 
large,  approziaately  2200,  in  comparison  to  the  other 
elements,  14.52  for  h  and  -6.4S  for  L.  Therefore, one  would 
expect  the  third  eigenvalue  to  be  very  sensitive  to  a  change 
in  the  thickness  and  less  sensitive  to  a  corresponding  change 
in  the  height  or  length.  Since  the  objective  function  is 
linear  with  respect  to  each  of  the  three  parameters,  a  large 
change  in  the  height  or  length  would  reduce  its  value  while 
having  less  effect  upon  the  eigenvalue  than  a  corresponding 
change  in  the  thickness.  Convergence  for  the  third 
eigenvalue  can  be  expected  to  a  much  lower  value  for  h  or  L 
with  t  remaining  approximately  the  same.  The  value  of  the 
objective  function  should  show  a  corresponding  decrease. 

CASE  1 

The  first  case  examined  cons t r a i ne d t h e  fundamental 
frequency  and  damping  coefficient  to  the  values  found  using 
the  nominal  thickness,  height,  and  length  given  by  Junkins  in 
[4]  and  using  the  gains  he  found  to  be  optimum.  The  first 
attempt  input  those  nominal  values  for  t,  h,  and  L  and  as 
expected,  very  little  change  occurred  as  the  program 
converged.  The  resulting  values  were  : 


.434780 

AM  = 

.000217 

.010411 

At  = 

.0000056 

.50 

Ah  - 

0 

4.0 

AL  - 

0 
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When  other  values  for  t,  h,  or  L  were  input,  the 
process  would  not  converge  and  a  feasible  solution  conld  not 
be  fonnd.  Since  the  inpnt  valnes  for  t,  h,  and  L  did  not 
return  valnes  for  g>j  and  gj  that  satisfied  the  constraints 
and  the  gradients  were  nearly  zero,  convergence  conld  not 
occur  and  after  10  iterations  a  feasible  solution  could  not 
be  found. 


CASE  2 


The  natural  frequency  and  damping  coefficient  of  the 
third  eigenvalue  were  constrained  to  be  equal  to  the  values 
found  using  the  nominal  t,  h,  and  L  and  the  gains  found  by 
Junkins.  Although  the  gradients  of  Xj  are  much  greater  than 
those  of  ,  the  program  still  would  not  reduce  the 
objective  function  any  significant  amount. 


iBJLSl 

t  =  .0104166 
h  =  0.50 
L  =  4.0 
M  =  .434997 

All  «  .00657 

Vhen  other  values  for  t,  h,  and  L  were  tried  the  program 
would  not  satisfy  the  constraints  and  a  feasible  solution 
could  not  be  found.  Different  values  were  tried  for  the 
constraint  thickness  parameter  and  larger  push  off  factors 


fiJLil 

.01044985 

.49152336 

3.9947516 

.428425 

1.5* 


(theta)  and  still  the  program  would  not  converge, 
example,  the  following  values  were  attempted 


As  an 


ill 


.5 


f 


t  -  .02  h  -  .5  L  -  4.0 

t  -  .015  h  -  .4  L  -  3.7 

Appendix  D  provides  s  sore  detailed  listing  of  soae  of  this 
ontpnt.  These  initial  values  do  not  satisfy  the  constraints 
and  represent  an  initial  non-feasible  solution,  which  CONMIN 
states  it  can  handle.  Valnes  for  CT  ■  -.1.  -.2,  -1.0.  -2.0 

t 

were  tried  with  the  above  values  for  t.  h,  and  L.  The 
constraint  thickness  specifies  a  region  for  which  the 
constraints  any  be  considered  satisfied,  so  by  increasing  the 
value  of  CT,  the  prograa  is  allowed  greater  flexibility  to 
solve  the  problea.  However,  even  with  the  abnoraally  large 
values  of  CT  (-1.0,  -2.0)  the  prograa  would  not  converge  to  a 
feasible  solution.  The  problea  was  due  to  CONMIN's  inability 
to  solve  equality  constraints  (although  the  user's  aanual 
states  that  it  can)  and  with  only  three  paraaeters  allowed  to 
change,  the  thickness  ,  height,  and  length,  there  was  not 
enough  freedoa  for  the  prograa  to  search  for  the  optiaua 
solution  and  satisfy  the  constraints. 

££££  1 

To  alleviate  the  probleas  experienced  with  the  equality 
constraint,  a  region  was  specified  in  which  to  constrain  the 
third  natural  frequency.  The  initial  value  of  was  4.38 
and  the  region  specified  by  the  new  constraints  was  ±  .2  froa 


coefficient  wee  removed  to  allow  greeter  freedom  of  movement 
for  the  program.  Therefore,  the  equality  constraint  imposed 
in  Case  2  was  replaced  by  an  inequality  constraint.  Several 
different  values  for  t,  h,  and  L  were  input  along  with 
different  values  for  the  constraint  thickness,  CT,  and 
DELFUN,  the  minimum  relative  change  in  the  objective  function 
to  indicate  convergence.  Different  values  for  the  push  off 
factor.  Theta,  were  also  combined  with  the  above.  Table  2 
lists  the  values  input  and  the  final  values  for  some  of  the 
cases  input.  All  converged  to  local  minimums  and  satisfied 
the  constraints  for  each  input.  A  more  detailed  listing  of 
these  results  are  presented  in  Appendix  E  .  Many  solutions 
were  found  which  satisfied  the  constraints  and  reduced  the 
mass  considerably.  This  was  due  to  the  steepness  of  the 
curve  representing  dX^/dt  and  also  the  second  partial 
derivative  of  M  was  equal  to  zero  for  all  three  parameters. 

Instability  in  the  fundamental  frequency  was 
encountered  when  the  constraint  thickness  parameter  was 
increased  to  -2.0.  This  increase  provided  additional 
flexibility  in  the  optimization  process  and  reduced  the  mass 
even  more  as  shown  in  Table  2.  However,  this  caused  Aj  to 
move  to  an  unstable  region  because  there  was  no  constraint 
imposed  to  keep  it  in  the  left  half  plane.  Since  more  than 
one  relative  minima  exists  and  some  of  these  represent 
unstable  solutions,  a  constraint  was  imposed  on  the  real  part 
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CASE  1 

Unstable  ▼ibrstion  st  the  fundaaental  frequency  is  not 
acceptable,  hence  another  constraint  was  iaposed  to  aaintain 
the  real  part  of  Xj  in  the  left  half  plane.  t»3  was  still 
constrained  to  reaain  in  the  previously  specified  region. 

Several  atteapts  were  aide  to  ainiaize  the  aass  with 
these  new  constraints.  Soae  converged,  but  did  not  reduce 
the  aass  aore  than  about  3%,  while  others  would  not  converge 
at  all.  Then  values  other  than  the  noainal  values  for  t,  h, 
and  L  were  input,  .0104166,  .5,  4.0  respectively,  the 
prograa  would  not  converge.  Several  atteapts  were  tried  with 
different  values  for  CT.  DELFUN,  DABFUN,  and  THETA,  but  the 
prograa  still  would  not  converge.  The  constraint  that  was 
not  satisfied  aost  often  was  the  constraint  on  «3.  The 
prograa  would  reduce  the  aass  and  aaintain  Xj  in  the  left 
half  plane,  but  «3  tended  towards  a  value  of  5. 0-5. 2. 
Therefore,  the  constraint  on  ej  was  changed  to  aaintain  it  in 
the  region  4.2  -  5.2  and  the  prograa  began  to  function  auch 
better . 

The  values  found  optiaua  in  Case  3  were  input  with  the 
new  constraints  and  the  results  are  shown  in  Table  3.  The 
only  convergence  problea  encountered  is  shown.  Several 
coabinations  of  values  for  CT,  DELFUN,  DABFUN,  and  THETA  were 
tried,  but  to  no  avail.  It  was  found  that  if  the  values  were 
not  close  to  a  feasible  solution  with  gradients  of  sufficient 
aagnitude  to  provide  direction,  the  prograa  would  satisfy  one 
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of  the  constraints  (maintain  o»3  in  the  specified  region),  bnt 
not  the  other  (maintain  in  the  left  half  plane).  The 
valves  input  for  t,  h,  and  L  represented  a  value  for  the  real 
part  of  X^  of  +1.815  which  evidently  was  too  far  off  to 
perait  convergence,  due  to  the  snail  gradients  of  X^. 

The  optiaua  values  were  found  to  be  as  shown  in  Table  3 
with  the  greatest  reduction  in  appendage  lass  of  60.7%.  A 
aore  detailed  listing  of  this  output  is  in  Appendix  F  along 
with  the  eigenvalues  and  eigenvectors  for  the  reconfigured 
structure. 


INPUT 

RESULT 

t 

h 

L 

CT 

t 

h 

L 

M 

X1 

.0104 

.5 

4  .0 

-.1 

.0127 

.494 

3.99 

.4238 

-.009 

.0104 

.5 

4.0 

-1 .0 

.0103 

.483 

3.99 

.4167 

-.010 

.012 

.1878 

3.58 

-.1 

Would  not 

satisfy 

+1.8 

.012 

.1878 

3.58 

-1 .0 

the 

constraint 

on 

+  1.8 

.012 

.1878 

3.58 

-2  .0 

X1 

+1 .8 

.0127 

.2695 

3.93 

-  .1 

.0125 

.2695 

3 .93 

.2772 

-.010 

.0127 

.2695 

3.93 

-1.0 

.0125 

.2695 

3.93 

.2772 

-.0095 

.012 

.3189 

3.94 

-.1 

.0118 

.319 

3.94 

.3108 

-.0099 

.0104 

.3 

3.8 

-1.0 

.0119 

.2997 

3.79 

.2841 

-.0074 

.0153 

.1434 

3.87 

-  .1 

.0156 

.1365 

3 .83 

.17015 

-.0089 

Table  3.  Results  for  Case  4. 


Further  study  should  be  done  to  expend  the  approach 


nsed  here  to  allow  for  varying  more  parameters.  The 
appendages  in  this  stndy  were  kept  structurally  identical 
which  reduced  the  number  of  degrees  of  freedom  in  the 
optimization  to  three  before  the  constraints  were  included. 
By  changing  the  appendages  to  pairs  of  structurally  identical 
beams,  the  parameter  vector  can  be  increased  by  a  factor  of 
two.  This  would  permit  greater  freedom  in  the  program  for 
optimization  of  the  parameters  and  minimization  of  the  total 
mass  while  constraining  the  frequencies.  The  cross-sectional 
areas  of  the  arms  were  kept  constant  from  the  root  to  the  tip 
in  this  report,  further  study  should  include  varying  the 
cross-sectional  area  such  that  tapered  beams  are  also 
included.  The  program  developed  here  to  call  CONMIN  should 
be  expanded  to  include  an  algorithm  to  optimize  the  closed 
loop  control  law  for  the  optimized  structure.  A  continuation 
procedure  could  be  used  as  an  iterative  technique  in  a 
similar  manner  to  that  used  by  Venkayya  in  [9]  or  by  Junkins 
in  [4].  The  program  should  begin  with  a  viable  solution  and 
iterate  until  the  mission  objectives  can  be  accomplished 
with  the  resulting  structure  and  control  law.  Another 


▼I.  CgBfifcBaiflM 


As  stated  by  Junkins  in  [4]  exact  pi aceaest  of 
eigenvalues  is  difficult  and  convergence  unlikely,  rather, 
specifying  a  particular  region  for  the  eigenvalue  leads  to  a 
■uch  better  algoritha  and  convergence.  By  reducing  the 
number  of  constraints  and  increasing  the  degrees  of  freedom 

in  the  problem,  the  program  converged  to  a  solution  which 

« 

minimized  the  mass  of  the  structure  while  keeping  the  first 
natural  frequency  in  a  stable  region  and  maintaining  the 
third  natural  frequency  in  a  specified  region.  By  changing 
the  constraint  from  an  equality  constraint  to  an  inequality 
constraint,  the  algorithm  functioned  much  better  and  more 
efficiently.  The  program  developed  in  this  report  can  be 
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APPENDIX  A 


EIGENVALUES  AND  EIGENVECTORS  Of  FREE  SYSTEM 

The  eigenvalues  and  eigenvectors  o£  the  free 
(uncontrolled)  system  were  found  by  using  program  PHICAL, 
which  is  listed  in  subroutine  form  in  Appendix  C.  to 
calculate  the  mass  and  stiffness  matrices.  A  modified 
version  of  program  MEIGLR,  listed  in  subroutine  form  in 
Appendix  C,  uses  the  mass  and  stiffness  matrices  to  calculate 
the  eigenvalues  and  corresponding  eigenvectors.  The  first  10 
eigenvalues  and  eigenvectors  are  listed.  The  first  one 
represents  the  rigid  body  mode,  while  the  remaining  nine 
represent  the  modes  published  by  Junkins  in  [4]  and  resemble 
those  of  a  cantilever  beam  with  a  lumped  mass  at  the  end.  A 
plot  of  these  modes  is  also  included.  The  actual  deflection 
in  the  v  direction  is  of  opposite  sign  to  computer  output. 


EIGENVALUE  1  IS 


. OOOOOOOOE+OO 
. OOOOOOOOE+OO 

ITS  FREQUENCY  IS  IN  RAD/ SEC 

. OOOOOOOOE+OO 
. OOOOOOOOE+OO 

AND  ITS  EIGENVECTOR  IS: 

. i OOOOOOOE+O 1  . OOOOOOOOE+OO 
. OOOOOOOOE+OO  . OOOOOOOOE+OO 
. OOOOOOOOE+OO  . OOOOOOOOE+OO 
. OOOOOOOOE+OO  . OOOOOOOOE+OO 
. OOOOOOOOE+OO  . OOOOOOOOE+OO 
. OOOOOOOOE+OO  . OOOOOOOOE+OO 
. OOOOOOOOE+OO  . OOOOOOOOE+OO 
. OOOOOOOOE+OO  . OOOOOOOOE+OO 
. OOOOOOOOE+OO  . OOOOOOOOE+OO 
. OOOOOOOOE+OO  . OOOOOOOOE+OO 
. OOOOOOOOE+OO  . OOOOOOOOE+OO 

EIGENVALUE  2  ISt 

. 191 23373E+02 
. OOOOOOOOE+OO 

ITS  FREQUENCY  IS  IN  RAD/SEC 

. 43730507E+01 
. OOOOOOOOE+OO 

AND  ITS  EIGENVECTOR  ISs 

- . 38933 111E-12  . OOOOOOOOE+OO 

- . 1 OOOOOOOE+O 1  . OOOOOOOOE+OO 

- . 29368652E-02  . OOOOOOOOE+OO 

- . 1 2880784E-0 1  . OOOOOOOOE+OO 

. 36673 1 67E-03  . OOOOOOOOE+OO 

-.1981 3906E-02  . OOOOOOOOE+OO 

. 1 OOOOOOOE+O 1  . OOOOOOOOE+OO 

. 29368652E-02  . OOOOOOOOE+OO 

. 1 2880784E-0 1  . OOOOOOOOE+OO 

36673 167E-03  .OOOOOOOOE+OO 

.1981 3906E-02  . OOOOOOOOE+OO 

EIGENVALUE  3  ISt 

. 62560428E+02 
. OOOOOOOOE+OO 
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ITS  FREQUENCY  IS  IN  RAD/SEC 


. 79095 130E+01 
. OOOOOOOOE+OO 

AND  ITS  EIGENVECTOR  ISj 

. 1 OOOOOOOE+O 1  . OOOOOOOOE+OO 

- . 93098297E+00  . OOOOOOOOE+OO 

.711301 04E-02  . OOOOOOOOE+OO 

1 1 74434 IE-01  . OOOOOOOOE+OO 

. 98375566E— 03  . OOOOOOOOE+OO 

18661096E-02  . OOOOOOOOE+OO 

- . 93098297E+00  . OOOOOOOOE+OO 

.711301 04E— 02  . OOOOOOOOE+OO 

- . 1 174434 1E-0 1  . OOOOOOOOE+OO 

. 98375566E-03  . OOOOOOOOE+OO 

- . 1 866 1 096E-02  . OOOOOOOOE+OO 

EIGENVALUE  4  ISt 

. 26752296E+04 
. OOOOOOOOE+OO 

ITS  FREQUENCY  IS  IN  RAD/SEC 

. 51722621E+02 
. OOOOOOOOE+OO 

AND  ITS  EIGENVECTOR  IS* 

- . 134 19288E-1 1  . OOOOOOOOE+OO 

- . 1 OOOOOOOE+O 1  . OOOOOOOOE+OO 

-.281 47448E+00  . OOOOOOOOE+OO 

.114791 28E-02  . OOOOOOOOE+OO 

- . 1 628366 1 E-0 1  . OOOOOOOOE+OO 

. 1 52407 07E— 02  . OOOOOOOOE+OO 

. 1 OOOOOOOE+O 1  . OOOOOOOOE+OO 

. 28 1 47448E+00  . OOOOOOOOE+OO 

-.11 479 l 28E-02  . OOOOOOOOE+OO 

. 1 628366 1 E-0 1  . OOOOOOOOE+OO 

- . 1 3240707E-02  . OOOOOOOOE+OO 

EIGENVALUE  5  I Si 

.281 77824E+04 
. OOOOOOOOE+OO 

ITS  FREQUENCY  IS  IN  RAD /SEC 

. 53082788E+02 
.OOOOOOOOE+OO 
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AND  ITS  EIGENVECTOR  ISi 


35977 156E-01 
.1 OOOOOOOE+O 1 
. 27 1929 1 OE+OO 
— . 20087475E-02 
. 15639824E-01 
1 5868837E— 02 
. 1 OOOOOOOE+O 1 
. 27 1 929 1 OE+OO 
20087 475E-02 
. 15639824E-01 
1 5868837E-02 

EIGENVALUE  6  IS: 

. 25682245E+05 
.  OOOOOOOOE+OO 


. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
.OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
.OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 


ITS  FREQUENCY  IS  IN  RAD/SEC 


. 1 602568 1E+03 
. OOOOOOOOE+OO 


AND  ITS  EIGENVECTOR  IS: 


•  267741 15E— 10 
. 1 OOOOOOOE+O 1 
— . 3025 1 849E+00 
- . 200 1 2248E+00 
. 997 65065E-02 
— . 22870335E— 0 1 
- . 99999999E+00 
. 3025 1 849E+00 
. 200 1 2248E+00 
- . 99765064E— 02 
. 22870335E-0 1 

EIGENVALUE  7  IS: 

. 25949087E+05 
. OOOOOOOOE+OO 

ITS  FREQUENCY  IS 

. 161 08720E+03 
. OOOOOOOOE+OO 


. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 
. OOOOOOOOE+OO 


IN  RAD/SEC 


AND  ITS  EIGENVECTOR  IS: 

. 85234335E-02  . OOOOOOOOE+OO 

. 1 OOOOOOOE+O 1  . OOOOOOOOE+OO 

-.31 23736 1 E+00  . OOOOOOOOE+OO 


- . 203406 75E+00  . OOOOOOOOE+OO 

. 10519143E-01  .OOOOOOOOE+OO 
- . 23206033E-0 1  . 00000000E+00 

. 99999999E+00  . OOOOOOOOE+OO 

3123736 1E+00  .  00000000E+00 

- . 20340673E+00  . OOOOOOOOE+OO 

. 10519143E-01  .OOOOOOOOE+OO 

- . 23206033E-0 l  . OOOOOOOOE+OO 

EIGENVALUE  8  IS: 

. 11410779E+06 
. OOOOOOOOE+OO 

% 

ITS  FREQUENCY  IS  IN  RAD/SEC 

. 33779845E+03 
. OOOOOOOOE+OO 

AND  ITS  EIGENVECTOR  IS: 

-. 1 68820 1 5E-09  . OOOOOOOOE+OO 

- . 99999996E+00  . OOOOOOOOE+OO 

. 73378994E+00  . OOOOOOOOE+OO 

-.41 778865E+00  . OOOOOOOOE+OO 

- . 43354098E+00  . OOOOOOOOE+OO 

. 52687S59E— O 1  . OOOOOOOOE+OO 

. 1 OOOOOOOE+O 1  . OOOOOOOOE+OO 

- . 73378997E+00  . OOOOOOOOE+OO 

. 4 1 778866E+00  . OOOOOOOOE+OO 

. 43334 1 OOE+OO  . OOOOOOOOE+OO 

5268736 IE-01  .OOOOOOOOE+OO 

EIGENVALUE  9  IS: 

. 11453340E+06 
. OOOOOOOOE+OO 

ITS  FREQUENCY  IS  IN  RAD/SEC 

. 33842783E+03 
. OOOOOOOOE+OO 

AND  ITS  EIGENVECTOR  IS: 

- . 82007604E-02  . OOOOOOOOE+OO 

. 99999996E+00  .OOOOOOOOE+OO 

- . 72262577E+00  . OOOOOOOOE+OO 

•41881917E+00  .OOOOOOOOE+OO 
. 43060762E+00  . OOOOOOOOE+OO 

- . 32829233E-0 1  . OOOOOOOOE+OO 

. 1 OOOOOOOE+O 1  . OOOOOOOOE+OO 

- . 72262380E+00  . OOOOOOOOE+OO 


. 4188191SE+00  . OOOOOOOOE+OO 

. 43060764E+00  . OOOOOOOOE+OO 

- „ 52829255E-0 1  . OOOOOOOOE+OO 

EIGENVALUE  10  ISt 

. 333676 18E+06 
. OOOOOOOOE+OO 

ITS  FREQUENCY  IS  IN  RAD/SEC 

. 377647 11E+03 
. OOOOOOOOE+OO 

AND  ITS  EIGENVECTOR  ISt 

.  133611 15E-09  .  OOOOOOOOE+OO 

- . 99999994E+00  . OOOOOOOOE+OO 

. 12120368E+00  . OOOOOOOOE+OO 

- . 4323472 1 E+00  . OOOOOOOOE+OO 

. 2797779 1 E+00  . OOOOOOOOE+OO 

.  4 1 077886E+00  .  OOOOOOOOE-KJO 

. 1 OOOOOOOE+O 1  . OOOOOOOOE+OO 

-.121 20369E+00  .  OOOOOOOOE+OO 

. 43234724E+00  . OOOOOOOOE+OO 

-.27977793E+00  .OOOOOOOOE+OO 
— . 41077889E+00  .OOOOOOOOE+OO 


y  yy*: 


Af?6Ki?II  1 


EIGENVALUES  AM  EJflMIM  Of  £ONTROLLEfi  SYSTEM 

The  eigenvalues  end  left  end  right  eigenvectors  were 
calculated  by  using  the  gains  fonnd  optimum  by  Jnnkina  in  [4] 
and  the  development  presented  in  Chapter  2  for  the  controlled 
system.  The  first  10  eigenvalues  and  corresponding 
eigenvectors  were  calculated  and  found  to  all  be  stable  (  in 
the  left  half  plane).  The  computer  program  used  to  perform 
the  calculations  is  described  in  Appendix  C  and  listed  in 
subroutine  form  as  MEIGLR.  The  optimum  control  lav  developed 
by  Junkins  resulted  in  low  frequencies  in  the  stable  left 
half  plane. 
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EIGENVALUE  1  IS: 


90228088E-0 1 
— . 31614593E+01 


AND  ITS  EIGENVECTOR  IS» 


,  68011813E— 12 
.  1 7276375E— 03 
.  70336592E-04 
.  17276372E-03 
,  70336580E-04 
.  38191026E—10 
,  10000000E+01 
,  35963626E+00 
,  lOOOOOOOE+Ol 
.35963626E+00 


-.67319834E-12 
. 1 8962959E— O 1 
. 68 1 96927E— 02 
1 8962939E-0 1 
-• 68196927E— 02 
-.36293787E-10 
. OOOOOOOOE+OO 
. 4337 1 307E-03 
. 18831664E-08 
— . 43371241E— 03 


ITS  LEFT  EIGENVECTOR  IS  : 


170681 15E+00 
. lOOOOOOOE+Ol 
. 26204397E+00 
. lOOOOOOOE+Ol 
. 26204397E+00 
■ . 56003427E— O 1 
. 12102 168E +00 
-.93076868E-03 
. 12102168E+00 
■.93076868E-03 


. 14464482E+00 
-.43882433E-12 
. 1 273 1 568E+00 
. OOOOOOOOE+OO 
. 12731568E+00 
-.321 18343E— 01 
. 60673223E-01 
1 686040 7E -03 
. 60673223E-0 1 
1 6860407E— 03 


EIGENVALUE  2  ISt 

— . 90228088E— O 1 
.3161 4593E+0 1 


AND  ITS  EIGENVECTOR  ISi 


. 68011813E— 12 
■.  17276373E-03 
70356392E-04 
. 1 727 6372E-03 
. 70356580E-04 
".38191026E— 10 
.10000000E+01 
.  33963626E+00 
•.  10000000E+01 
■ .  33963626E+00 


. 673 1 9834E— 1 2 
-. 18962959E-0 1 
68196927E-02 
. 18962939E— 0 1 
. 68196927E-02 
. 36293787E-10 
. OOOOOOOOE+OO 
-. 43371307E-03 
-. 18831 664E-08 
.43371 24 IE-03 


ITS  LEFT  EIGENVECTOR  IS  * 


-. 170681 15E+00 
. 1 OOOOOOOE+O 1 
• 26204397E+00 
. 1 OOOOOOOE+O 1 
. 26204397E+00 
- . 56003427E-0 1 
. 1 21021 68E+00 
-.95076868E-03 
. 12102168E+00 
- . 95076868E-03 


1 44644S2E+00 
. 4388243SE— 1 2 
1 273 1 368E+00 
. OOOOOOOOE+OO 
12731568E+00 
. 32 1 1 8343E-0 1 
— . 60673223E— 01 
. 1 686040 7E— 03 
-.60673225E-01 
. 1 686040 7E— 03 


EIGENVALUE  3  ISt 

1 3384354E+00 
— . 4387 6209E+0 1 


AND  ITS  EIGENVECTOR  ISt 


. 40064767E— 1 1 
69461 113E-02 
47790883E-03 
. 69461 113E-02 
. 47790883E— 03 
19078769E— 10 
. 1 OOOOOOOE+O 1 
13177105E— Ol 
1 OOOOOOOE+O 1 
. 13177103E— 01 


- . 37278 1 03E~ 1 1 
. 22770205E+00 
— . 30 1 73248E-02 
22770203E+00 
.301 78248E-02 
1 6823943E- 1 0 
. 10981794E—14 
. 23008033E—02 
•  32901435E'-1 1 
— . 23008033E—02 


ITS  LEFT  EIGENVECTOR  IS  l 


— • 840 1 0789E—02 
94582793E+00 
-.26137172E+00 
. 1 OOOOOOOE+O 1 
. 27293433E+00 
— . 22747278E— 02 
87306233E-01 
. 1 383232 IE-02 
. 9339643 1E-0 l 
1 453 1 629E-02 


.911 16440E—02 
-.4731 937 8E-02 
1 4350268E+00 
. OOOOOOOOE+OO 
. 15131088E+00 
-. 1 60 1 2045E-02 
— . 46890308E-0 1 
. 1 8880546E-03 
.511 1662SE— 01 
— . 21337217E-03 


EIGENVALUE  4  ISt 


13384354E+00 

43876209E+01 


AND  ITS  EIGENVECTOR  IS* 


.40064767E-11 
-.69461 11 3E-02 
47790885E— 03 
. 6946 1 1 13E— 02 
. 47790885E-03 
19078769E— 10 
. 1 OOOOOOOE+O 1 
— . 13177105E— 01 
1 OOOOOOOE+O 1 
. 13177105E— 01 


. 37278 105E- 11 
-.22770 205E+O0 
. 30 1 7S248E— 02 
.22770205E+00 
-. 30 1 78248E-02 
. 1 6823943E— 1 O 
-. 10981794E— 14 
- . 2S0080S3E— 02 
— . 32901455E— 1 1 
. 2S008053E-02 


ITS  LEFT  EIGENVECTOR  IS  * 


■  . 840 1 0789E— 02 
94582793E+00 
■.26137172E+00 
. 1 OOOOOOOE+O 1 
. 27293433E+00 
■.  22747278E— 02 
>.  87306233E-01 
. 1 383232 IE-02 
.93396431E-01 
•.  14351629E-02 


— .91116440E— 02 
. 473 1 9578E-02 
. 1 4330268E+00 
. OOOOOOOOE+OO 
—  .  13131088E+00 
. 1 60 1 2045E-02 
. 46890308E-01 
-. 1 8880546E— 03 
-.51 1 16625E— 01 
•21357217E— 03 


EIGENVALUE  3  IS* 

-. 15949736E+00 
— . 4756 1360E+0 1 


AND  ITS  EIGENVECTOR  IS* 


- . 9020 1 333E— 02 
-. 21421430E— 02 
- . 24 1 83730E— 04 
-. 20928402E-02 
— . 23036809E— 04 
. 1 OOOOOOOE+O 1 
-.34031291E+00 
-. 1 4736053E— 0 1 
-. 34280870E+00 
-.  1 4769564E-0 1 


. 316052 19E+00 
-. 1 0770338E+00 
— . 4668 171 6E— 02 
-. 1 0849342E+00 
-. 46724 133E-02 
- . 1 4370798E- 1 4 
. 1 6490349E-0 1 
. 49763608E-03 
. 1 6405383E-0 1 
. 49447608E-03 


ITS  LEFT  EIGENVECTOR  IS  * 


. 39933862E+00 
. 1 OOOOOOOE+O 1 
. 17964978E+00 
. 1 OOOOOOOE+O 1 
. 17964978E+00 


. 3483307 1E+00 
— . 48223966E— 12 
. 1 6696430E+00 
. 17729399E-14 
. 1 6696430E+00 
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- . 805 12903E-0 1  - . 755276 1 7E-0 1 

. 1024621 1E+00  . 93405672E-0 1 

- . l 3506533E-02  -.71 085S98E-03 

. 1024621 lE+OO  . 93405672E— O 1 

-. 13S06535E-02  -. 71085598E-03 


EIGENVALUE  6  IS* 

15949736E+00 
. 47S6 1 360E+0 1 


AND  ITS  EIGENVECTOR  IS* 

9020 1 333E-02  -. 31605219E+00 

21421430E-02  . 10770538E+00 

— . 24 1 83730E— 04  . 46681716E-02 

- . 20928402E-02  . 1 0849342E+00 

- . 23056809E— 04  . 46724 1 33E-02 

. 1 OOOOOOOE+O 1  . 14370798E-14 

-.34031 29 1E+00  1 6490349E-0 1 

- . 1 4736053E— O 1  - . 49765608E-03 

— . 34280870E+00  -. 16405S83E-01 

-. 14769564E-01  49447608E-03 

ITS  LEFT  EIGENVECTOR  IS  s 
-.39933862E+00  3483307 1E+00 

. 1 OOOOOOOE+O 1  . 48223966E— 1 2 

. 1 7964978E+00  - . 1 6696430E+00 

. 1 OOOOOOOE+O 1  1 7729399E— 1 4 

. 1 7964978E+00  - . 1 6696430E+00 

80512903E— 01  . 75327617E-01 

. 1024621 lE+OO  -. 93403672E-01 
- . 1 3506533E-02  . 7 1 085598E-03 

. 1024621 1E+00  -. 93403672E-0 1 

- . 1 3506535E— 02  . 7 1 083598E-03 

EIGENVALUE  7  IS: 

48040 152E+00 
. 327300 10E+02 


AND  ITS  EIGENVECTOR  IS* 

-. 7042986 IE-02  -. 21001832E+00 

— .61 13S169E-03  . 13157677E+00 

. 17133024E-03  . 39223396E-02 

-. 7670 1 797E-03  . 13024501E+00 

. 1 70033 1 8E-03  . 3950 1 266E-02 

. 10000000E+01  . 30518120E— 14 

- . 62569949E+00  - . 2389396 1 E-0 1 


1868231 IE-01  . 1 8936222E-03 

-.61934063E+00  24421 776E-01 

18814462E-01  . 17876 127E-03 

ITS  LEFT  EIGENVECTOR  IS  : 

-. 15044437E-06  -.43499205E-07 

- .  99988239E+00  .  28759423E-05 

- . 35959628E+00  . 23707392E-05 

. iOOOOOOOE+O 1  - . 80986904E- 1 3 

. 33963852E+00  -. 15594396E-05 

-.271 68420E-08  . 279 1 8 1 20E-08 

1 83933S6E-03  . 18958907E-01 

— . 61921343E— 04  . 68 1 83756E-02 

. 1 8403 1 OOE— 03  18961 137E-01 

. 61947 106E— 04  -. 68191764E-02 


EIGENVALUE  8  IS: 

48040 132E+00 
527300 10E+02 


AND  ITS  EIGENVECTOR  IS: 

7042986 1 E— 02  . 21001832E+00 

61 138169E-03  1 31376 77E+00 

.171 33024E— 03  - . 39223396E-02 

— . 76701797E— 03  -. 1 302430 lE+OO 

. 1 70033 1 8E— 03  39301266E-02 

. lOOOOOOOE+Ol  30318120E-14 

- . 62369949E+00  . 2389396 1 E-0 1 

1868251 IE-01  -. 1 8936222E-03 

- . 6 1934063E+00  . 24421776E-01 

18814462E-01  -. 17876 127E-03 

ITS  LEFT  EIGENVECTOR  IS  : 

- . 1 5044437E— 06  . 43499205E-07 

- . 99988239E+00  - . 28759423E-05 

- . 33939628E+00  - . 25707392E-05 

. 10000000E+0 1  . 80986904E- 1 5 

. 33963832E+00  . 15594396E-03 

- . 27 1 68420E— 08  -. 27918120E-08 

- . 1 8395336E-03  - . 1 8958907E-0 1 

61921343E-04  - . 68 1 83736E-02 

. 18403 100E-03  . 18961 137E-01 

.61947106E— 04  . 68191764E-02 


EIGENVALUE  9  IS: 

- . 49446867E+00 
. 539 1 7505E+02 


AND  ITS  EIGENVECTOR  IS* 


13121 624E-04 
17007586E-03 
-.71610986E-04 
17010004E-03 
— . 71618821E— 04 
-.39612890E-01 
. 10000000E+01 
.34911367E+00 
. 99989 129E+00 
. 34907 468E+00 


.  7348 1 477E—03 
18545293E-01 
— .647 430 29 E- 02 
18543276E— O 1 
64735798E-02 
1070828 IE-02 
. 69625827E— 1 5 
-.65974577E-03 
23008404E-05 
6605257 4E-03 


ITS  LEFT  EIGENVECTOR  IS  : 


- . 36828520E-0 1 
. 1 OOOOOOOE+O 1 
•34997737E+00 
. 1 OOOOOOOE+O 1 
. 34997737E+00 
270429 12E-04 
. 1 9999 138E— 03 
. 591 58834E— 04 
. 19999 138E-03 
. 59 1 58834E— 04 


-.72612925E-05 
1 0866528E- 1 1 
— . 37209655E— 05 
. OOOOOOOOE+OO 
-. 3720965 IE-05 
. 68269453E-03 
1 8542003E—0 1 
- . 64893800E— 02 
1 8542003E— 0 1 
64893800E-02 


EIGENVALUE  10  IS* 

49446867E+00 
- . 539 1 7505E+02 


AND  ITS  EIGENVECTOR  IS: 


13121624E-04 
17007586E-03 
— . 71610986E— 04 
1701 0004E— 03 
71618821E-04 
— . 39612890E-01 
. 1 OOOOOOOE+O 1 
. 3491 1367E+00 
. 99989 129E+00 
. 34907 468E+00 


73481477E-03 
. 1 8545293E-0 1 
. 64743029E-02 
. 1 854327 6E— O 1 
. 64735798E-02 
. 1 070828 IE-02 
— . 69625827E- 1 5 
. 6597 4577E-03 
. 23008404E— 05 
• 66052574E— 03 


ITS  LEFT  EIGENVECTOR  IS  s 
- . 36828520E—0 1  .7261 2925E-OS 

. 1 OOOOOOOE+O 1  . 1 0866528E- 1 1 

. 34997737E+00  . 37209655E-05 

. 1 OOOOOOOE+O 1  . OOOOOOOOE+OO 

. 34997737E+00  . 3720965 1 E-05 

270429 12E-04  68269453E-03 

. 1 9999 1 38E-03  . 1 8542003E-0 1 

. 59 1 58834E— 04  . 64893800E-02 

. 1 9999 1 38E-03  . 1 8542003E-0 1 

. 59 1 58834E-04  . 64893800E-02 


APPENDIX  C 


COMPUTER  PROGRAMS 

The  following  computer  program  contains  the  programs 
used  in  this  report  in  subroutine  form, except  for  subroutine 
CONMIN  which  can  be  obtained  through  NASA  or  Dr.Khot,  VPAFB. 

Program  MASMIN  initializes  all  parameters  required  by 
CONMIN  and  dimensions  the  arrays  common  to  CONMIN.  It 
prompts  the  user  for  values  for  DELFUN,  DABFUN,  CT,  and  THETA 
and  then  asks  for  the  initial  values  for  t,  h,  and  L.  The 
lower  and  upper  bounds  on  the  decision  variables  are  set  and 
CONMIN  called  to  begin  the  process.  Subroutine  CALOBJ  is 
called  to  calculate  the  objective  function,  gradients  of  the 
objective  function,  constraints  and  gradients  of  the 
constraints.  File  MINMAS  is  passed  the  final  values  of  the 
objective  function  and  decision  variables. 

Subroutine  CALOBJ  calculates  the  gradients  of  the 
constraints,  which  were  w  and  §  of  a  specified  eigenvalue,  by 
calling  subroutine  PHICAL.  Matrix  A,  which  is  common  to 
CONMIN,  contains  the  gradients  of  the  constraints  while  array 
DF  contains  the  gradients  of  the  objective  function. 

Subroutine  PHICAL  calculates  the  mass  and  stiffness 
matrices  for  specified  values  of  t,  h,  and  L.  This  data  is 
stored  in  files  'MMASS'  and  'MSTIFF'  for  separate  access. 
PHICAL  calls  subroutine  CALC1  which  computes  values  for  $(z), 
+  '(z),  integral  of  ^  fj,  integral  of  x^j^z),  and  the 
integral  of  t'it’i  numerically.  These  values  may  be  printed 


s 


» 


p 


i 


I 


oat  by  the  user  if  desired  end  are  used  to  form  the  mass  and 
stiffness  matrices. 

Subroutine  NEIGLR  is  passed  the  mass  and  stiffness 
matrices,  the  length  of  the  appendages,  and  the  number  of 
modes  being  examined  by  PHICAL  and  it  uses  this  data  to  form 
matrices  A  and  B.  An  IMSL  subroutine  called  EIGZF  is  passed 
matrices  A  and  B  which  solves  for  the  eigenvalues  and 
eigenvectors.  The  transpose  of  A  and  B.  are  also  passed  to 
EIGZF  to  return  the  left  eigenvectors.  MEIGLR  forms  the  A 
and  B  matrices  by  calculating  the  B(called  BZ  in  the  program) 
and  H  matrices  as  shown  in  Chapter  2  and  calling  subroutine 
Gain  which  contains  the  gain  valnes.  The  matrices  are 
assembled  with  the  previously  calculated  mass  and  stiffness 
matrices.  Matrix  A  and  B  may  also  be  printed  if  desired 
along  with  the  eigenvalues  and  eigenvectors. 

Subroutine  CALCBA  is  passed  the  specified  or  input 
values  for  t,  h,  L,  the  number  of  modes  being  examined,  the 
values  of  the  integrals  calculated  by  PHICAL  along  with  the 
eigenvalues  and  eigenvectors.  It  outputs  the  gradients  of  <■> 
and  £  for  the  specified  eigenvalue.  In  other  words  if  the 
third  eigenvalue  is  being  used  as  a  constraint,  then  the  user 
should  output  the  gradients  of  and  (3.  The  gradients  of 
the  mass  and  stiffness  matrices  along  with  Bz  and  H  are 
computed  for  each  structural  parameter  and  used  to  form  the 
gradients  of  the  A  and  B  matrices. 


The  gradients  of  A  and  B  are  passed  to  subroutine 


GRADNT  which  forms  the  gradients  of  the  eigenvalues.  This  is 
accomplished  in  accordance  with  the  development  in  Chapter  3 
using  the  partial  derivatives  of  A  and  B  (Chapter  3  refers  to 
matrix  B  as  D  to  avoid  confusion  with  the  control  matrix  B). 
Once  the  gradients  of  the  eigenvalues  are  found,  the 
gradients  of  <o  and  $  are  easily  formed  and  are  referred  to  as 
PW3  and  PZ3.  These  represent  3  element  arrays  containing  the 
partial  derivative  of  <1)3  w.r.t.  t,  h,  and  L  and  the  partial 
derivative  of  £3  w.r.t.  t,  h,  and  L,  respectively. 

This  information  is  passed  back  to  subroutine  CALOBJ 
which  aids  CONMIN  in  its  constrained  minimization  problem. 
The  gradients  of  U3  and  £3  along  with  those  of  the  objective 
function  provide  direction  for  the  method  of  Feasible 
Directions  used  by  CONMIN. 


PROGRAM  M ASM IN 

COMMON/CNMN i /DELFUN, DABFUN, FDCH, FDCHM, CT, CTMIN, CTL, CTLMIN, 

1 ALPHA X , ABOBJ 1 , THETA, OBJ, NDV, NOON, NS IDE, IPRINT, NFDG, 

2NSCAL , L I NOB J , ITMAX, ITRM, ICNDIR, IGOTO,NAC, INFO, INF06, ITER 
COMMQN/CONSAV/REAL (50) , I NT (25) 

DIMENSION  X  (50) ,  VLB  (50)  ,  VUB  (50)  ,6(1 100) ,  SCAL  (50) ,  DF  (50) 

1A (5,51) ,S(51) ,GI (1100) ,G2(I100) ,B (51,51) ,C(51) , 

2 I SC ( 1000) , IC(51) ,MS1 (102) 

C 

C  DELFUN  AND  DABFUN  SPECIFY  THE  TOLERANCE  ON  THE  OBJECTIVE 

C  FUNCTION,  CT  IS  THE  CONSTRAINT  THICKNESS  PARAMETER  AND 
C  SPECIFIES  THE  TOLERANCE  ON  SATISFYING  THE  CONSTRAINTS.  THETA 
C  IS  THE  PUSH  OFF  FACTOR  AND  IPRINT  ALLOWS  THE  USER  TO  SPECIFY 
C  HOW  MUCH  OUTPUT  HE  WANTS  FROM  CONMIN. 


PRINT*, * PLEASE  INPUT  DELFUN, DABFUN, CT, THETA, IPRINT* 
READ*, DELFUN, DABFUN, CT, THETA, IPRINT 
PRINT*, 'PLEASE  INPUT  THE  VALUES  FOR  T,H,L* 

DO  15  1*1 , 3 
READ* , X ( I ) 

CONTINUE 

FDCH  »  0.0 
FDCHM  *0.0 
CTMIN  *  0.0 
CTL  *  0.0 
CTLMIN  -  0.0 
ALPHA X  *  0.1 
ABOBJ 1  *0.1 
NDV  -  3 
NCON  *  2  . 

NFDG  »  1 
NSCAL  *  O 
LINOBJ  -  O 
ITMAX  *  40 
ITRM  *  O 
ICNDIR  -  0 
NS IDE  «  6 
I GOTO  -  O 
N1  »  5 
N2  *  S 
N3  *  51 
N4  *  31 
N3  »  102 


DO  10  I  *  1 , NCON 


ISC(I)  -  0 
10  CONTINUE 

DO  20  I  -  1 ,NDV 
VLB(I)-. 000001 
VUB ( I ) »1 . 0E+ 1 0 


77 


U  U  UwUUOUU 


20  CONTINUE 


OPEN (UNIT-9,  FILE-* MINMAS’ , ACCESS-’ SEQUENTIAL* , STATUS- ’ NEW * ) 
REWIND (UNIT-9) 

50  CALL  CONMIN  <X,  VLB,  VUB,  G,  SCAL, DF, A,  S,  61 , 62,  B,  C,  ISC, 

1  IC, MSI , Nl, N2, N3, N4, N5) 

IF  (IGOTO  .EQ.  0)  GO  TO  200 
IF  (INFO  .EQ.  2)  GO  TO  150 
CALL  CALOBJ (X,G,DF, A, IC, N3) 

GO  TO  50 

150  CALL  CALOBJ <X,G,DF, A, IC,N3> 

GO  TO  50 

200  CALL  CALOBJ <X, G,DF, A, IC,N3> 

WRITE (9, 185)  X ( 1 ) ,X<2> ,  X  (3) 

185  FORMAT < 5X , 7HX  < 1 )  »  , E15. 8, 5X , 7HX (2)  »  ,E15.8,5X, 

17HX (3)  *  , E15.8/) 

WRITE (9, 187) OBJ, G( 1) ,6(2) 

187  FORMAT <5X,21H0BJECTIVE  FUNCTION  «  ,E15.8,/5X, 

121 HCONSTRA I NT  1  -  ,E15.8,/5X, 

22 1 HCONSTR A I NT  2  -  ,E15.8/) 

ENDFILE (UNIT— 9) 

1000  STOP 

END 

************************************************************* 

************************************************************* 
SUBROUTINE  CALOBJ <X,G,DF, A, IC,N3) 

COMMON / CNMN 1 /DELFUN, DABFUN, FDCH, FDCHM, CT, CTMIN, CTL, CTLMIN, 

1 ALPHA X , ABOBJ 1 , THETA, OBJ , NDV, NCON, NSIDE, IPRINT, NFDG, 

2NSCAL, LINOBJ , ITMAX, ITRM, ICNDIR, IGOTO, NAC, INFO, INFOG, ITER 
DIMENSION  X (50) , G ( 1 100) , DF (50) , IC (51 ) , A (5, 51 ) 

INTEGER  I 
COMPLEX  LBDA(IO) 

REAL  MD,T,H,L,PW3(3) ,PZ3(3) ,W3,Z3,Z 

**************************************************** 


THIS  CALCULATES  THE  OBJECTIVE  FUNCTION  AND  GRADIENTS 


200  FORMAT ( 5X , * T  -  * , E15. 8, 5X , * H  «  ’ , E15. 8, 5X, * L  -  *,E15.S/) 

C 

T-X  ( 1 ) 

H«X(2> 

L-X(3) 

C 

MD-5.22 

CALL  PHICAL  (.T ,  H,  L,  PW3,  PZ3,  LBDA) 

C  INPUT  IS  T,H,L 


IF  (INFQ.EQ.2)G0  TO  10 


C 

C  CALCULATE  OBJECTIVE  FUNCTION  AND  GRADIENTS 
C 

0BJ>4*MD*T*H*L 

WRITE(#, *> 'THIS  IS  ITERATION  » , ITER 
WRITE (*, *) *  THIS  IS  THE  OBJECTIVE  FUNCTION  * ,OBJ 
WRITE (*,*> ’THIS  IS  LAMBDA  ’,LBDA(3> 

WRITE(*,200>  X(l> ,X(2> ,X(3> 

C 

C  CALCULATE  GRADIENTS  6(1)  AND  G(2> 

C  M3  IS  RE-COMPUTED  VALUE  OF  OMEGA 

C  Z3  IS  RE-COMPUTED  VALUE  OF  ZETA,  FOR  THIRD  EIGENVALUE 

C 

W3  “SORT ( (REAL (LBDA (3> ) ) **2  +  ( AIMAG (LBDA (3) ) ) **2) 
Z*REAL (LBDA (3) ) 

Z3— Z/W3 
G ( 1 ) *W3— 4 . 6000 
G(2)*  4.2-W3 
C 

IF  (INFO.EQ. 1)  RETURN 
10  CONTINUE 

C 

DF ( 1 ) »4*MD*H*L 
DF (2) «4*MD*T*L 
DF(3)-4*MD*T*H 
C 

NAC*0 

IF  (G(l).LT.CT)  GO  TO  20 
NAC»1 
IC(1)»1 

A(i, 1)»-PW3(1) 

A  (2,  1>—  PW3(2) 

A (3, i)»— PW3(3) 

C 

20  IF  (G(2) .LT.CT)  RETURN 

NAC-NAC+1 

IF  (NAC.EQ.N3)  RETURN 
IC(NAC)»2 
A ( 1 , NAC) »PW3 ( 1 ) 

A(2,NAC)-PW3(2) 

A(3,NAC)»PW3(3) 

C 

C 

END 

C 

*************************************************************** 


SUBROUT I NE  PHI CAL (T ,  H,  L, PW3 , P Z3 , LBOA > 

INTEGER  I,J,N,K 

REAL  IX,PW3(3) ,PZ3(3) 

REAL  PHI <10) ,PHISQ<10, 10)  , IPHISQ(10, 10) , PHI1 (10) , PHI 1SQ (10, 10) 
REAL  IXPHI ( 10) , TEMP1A (10, 10) , TEMP1B ( 10, 10) , TEMP2A ( 10) , TEMP2B ( 10) 
REAL  M(21,21) ,C1,C, I2C, 1 1 ,M2,D,  T,  H,  AX,  E 
REAL  HZ <5, 5) 

COMPLEX  LIGVCT ( 10, 10) , RIGVCT ( 10, 10) ,LBDA<10) 

REAL  STIFF (21, 21) 

C  IS  THE  MOMENT  OF  INERTIA  OF  THE  HUB 

I2C  IS  THE  ROTARY  MOMENT  OF  INERTIA  OF  THE  PT  MASS 

D  IS  THE  DENSITY  OF  THE  ARMS 

II  IS  THE  TRANSLATED  MOMENT  OF  INERTIA  OF  THE  ARMS 
T  IS  THE  THICKNESS  OF  THE  ARMS 
H  IS  THE  HEIGHT  OF  THE  ARMS 
L  IS  THE  LENGTH  OF  THE  ARMS 
E  IS  THE  MODULUS  OF  ELASTICITY  OF  THE  ARMS 
IX  IS  THE  AREA  MOMENT  OF  INERTIA  OF  THE  ARMS 
REAL  R,  X,L, IPHI2< 10, 10) , TEMP3A (10, 10) , TEMP3B ( 10, 10) 

N-2 

R«1 

D*5.22 
M2*. 156941 
I2C-.0018 
E*l.584E9 
X-R 

CALL  CALC1 (X,L, N, PHI , PHISQ, TEMPI A, PHI 1 ,PHI 1SQ, TEMP2A, TEMP3A) 
X-R+L 

CALL  CALC1 (X,L, N, PHI, PHISQ, TEMP1B, PHI 1 , PHI ISO, TEMP2B, TEMP3B) 

DO  40  1*1, N 
DO  50  J»1,N 

I PHISQ ( I , J ) -TEMP1B (I, J) -TEMPI A (I ,  J) 

I XPH  Id)  -TEMP2B  ( I )  -TEMP2A  ( I ) 

IPHI2 ( I , J) »TEMP3B (I, J) -TEMP3A ( I , J ) 

0  CONTINUE 
CONTINUE 

TO  PRINT  THE  PHI , PHISQ, IPHISQ  ETC.  MATRICES  THEN  REMOVE  THE 
IF  STATEMENT. 

IF ( I . EQ. Ill) THEN 

PRINT*, ’THE  FOLLOWING  IS  PHI’ 

PRINT  100, (PHI (I), 1-1, N) 

PRINT*, ’THE  FOLLOWING  IS  PHISQ’ 

DO  60  1*1, N 

PRINT  100,  (PHISQd,  J)  ,  J*1,N) 

0  CONTINUE 

PRINT*, ’THIS  IS  IPHISQ’ 

DO  70  1*1, N 

PRINT  100,  (IPHISQd, J)  , J*t,N) 

0  CONTINUE 

PRINT*, ’THIS  IS  PHI l’ 


PRINT  100, (PHIl (I) , I-1,N> 

PRINT*, 'THIS  IS  PHI 1SQ' 

DO  80  I-1,N 

PRINT  100, (PHI1SQ(I, J) , J-1,N) 

80  CONTINUE 

PRINT*,* THIS  IS  IXPHI * 

PRINT  100, (IXPHI (I) , I-l,N> 

PRINT*, ’THE  FOLLOWING  RESULTS  ARE  IPHI2* 

PRINT*, *  * 

DO  85  1*1 ,  N 

PRINT  100, (IPHI2(I, J) ,J-1,N> 

85  CONTINUE 
ENDIF 

C-8.0 

AX»D*T*H 

I 1* (AX/3) * ( (R+L) **3-R**3> 

Cl* (C/2) +2*1 1+2*M2*  <R+L) **2 
IX*(1/12.0) *H*T**3 
ti  ( 1 , 1 )  *2*  (C1  +  I2C) 

STIFF ( 1, 1)*0.0 

DO  110  1*1, N 

M(l, I+l) *2* (AX* IXPHI ( I ) + (R+L) *M2*PHI ( I ) +I2C*PHI 1 (I) ) 

M ( 1 , 1+N+l >  — M ( 1 , 1+1 ) 

M ( I+l , 1) — M (1,1+1) 

M(I+N+1,1)»M(1,I+1) 

STIFF (1, I+l) -0.0 
STIFF (1,I+N+1) *0.0 
STIFF (1+1,1) -0.0 
STIFF ( I+N+l, 1) -0.0 
110  CONTINUE 

DO  120  1*1, N 
DO  130  J— 1 , N 

M(I+1,J+1) *2* (AX* I PH I SO ( I , J ) +M2*PHISQ ( I , J) +I2C*PHI ISO ( I , J ) ) 
M ( I+N+l , J+N+l ) -M ( I+l , J+l ) 

M(I+1,  J+N+D-0 
lid+N+t,  J+l)— 0 

STIFF(I+1, J+1)«2*(E*IX*IPHI2(I, J) ) 

STIFF (I+l+N, J+l) -0.0 
STIFF ( I+l , J+l+N) -0. 0 

STIFF ( I+l+N, J+l+N) -2* (E* IX*IPHI2 ( I , J) ) 

130  CONTINUE 

120  CONTINUE 

100  FORMAT (8 (E16.  8) ) 

CALL  ME I GLR ( L , N , M, ST I FF , L I GVCT , R I GVCT , LBDA ) 

CALL  CALCBA(T,H,L,N, IXPHI, IPHISQ, IPHI2, LIGVCT, 

1  RIGVCT, LBDA, PW3, PZ3) 

END 
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SUBROUT I NE  CALC 1  ( Z ,  L ,  N , TEMP 1 , TEMP2 , MP , TEMP4 , TEMPS , TEMP 6 
REAL  X,L, TEMPI *TBW?7EMP2 ( 10, 10) ,  MP < 10, 10) , TEMP4 < 10) 
REAL  MP1,MP2,MP3,MP4,MP5, Z, R 
INTEGER  I , J, N 

REAL  TEMPS ( 10, 10) , TEMP6( 10) ,  A,  B,  C,  D,  F,  TEMP7 ( 10, 10) 

R-1.0 

X-Z-R 

PI-3. 141593 
DO  10  I»1,N 
DO  20  J»1,N 
A-I*PI/L 
B-J*PI/L 
C»(-i)**(I+l) 

D«(-1)**(J+1) 

IF  ( I .NE. J)  THEN 

F- ( 1 / (2$ ( A-B) > ) *SIN ( (A-B) #X ) * ( 1 / (2* ( A+B) ) ) 

L  *SIN< (A+B) *X) 

ELSE 

F-.5*X  +  ( 1/  (4*A)  )  *SIN <2*AtX") 

ENDIF 

TEMPI  ( I )  -1-COS  <  A*X )  -*•.  S*C*  <  A*X )  **2 

TEMP2 ( I , J ) -1-COS  <B*X ) +. 5*D*  (B*X ) **2  -COS (A*X) +COS (A*X  > 
L  *COS  <B*X )  — 

L  . 5*D* (B*X) **2*C0S ( A*X ) +  . 5*C* <A*X) **2-. 5*C* ( A*X ) **2 

L  *COS <B*X) 

L  +. 25*C*D* ( A*B) **2*X**4 

MP1-X-(1/B)*SIN(B*X)+(1/6.0)*D*B*B*X**3-(1/A) 

L  *SIN(A*X) 

MP2— .5*D*B*B* 

L  (  (1/A**2)*2*X*C0S(A*X)+(1/A**3)*(  <A*X)**2-2) 

L  *SIN ( A*X) ) 

MP3— (1/6.0) *C*A*A*X**3 

MP4— .3*C*A*A*( (1/B**2) *2*X*C0S (B*X) + ( 1/B**3) 

L  *( (BtX)**2-2)*SIN(B*X) ) 

MP3- ( 1/20.0) *C*D*A*A*B*B*X**5  +  F 
MP ( I , J ) -MP1+MP2+MP3+MP4+MP5 
TEMP4 (I) -A*SIN (A*X) +C*A*A*X 

TEMPS (I, J)»A*B*SIN(A*X)*SIN(B*X)+A*D*B*B*X*SIN(A*X) 

L  +A*A*B*C*X*SIN(B*X)+A*A*B*B*C*D*X*X 

TEMP6(I)».3*Z*Z-(1/A**2) *COS(A*X)-(Z/A) *SIN(A*X) 

L  +.3*C* ( A**2) * (.25* (Z**4> -(2/3.0) *R*Z**3+.5*R*R*Z*Z) 
TEMP7(I,  J)»(  (A*B)**2)*(F-MC/B)*SIN(B*X)  +  (D/A)*SIN(A*X) 
L  +C*D*X) 

CONTINUE 

CONTINUE 


SUBROUTINE  MEIGLR (L, P, M, S, LIGVCT, Z^EIGENV) 
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******************************************************************* 

******************************************************************* 

THIS  PROGRAM  SOLVES  FOR  THE  EIGENVALUES  AND  EIGENVECTORS  BY 
SOLVING  THE  PROBLEM  AX  «  LAMBDA  BX.  A  IS  FORMED  FROM  THE 
STIFFNESS  AND  THE  CONTROL  MATRICES,  WHILE  B  IS  FORMED  FROM 
THE  MASS  MATRIX. 

Z  CONTAINS  THE  RIGHT  EIGENVECTOR  ,  LIGVCT  CONTAINS  THE  LEFT 
EIGENVECTOR,  AND  EIGENV  IS  THE  EIGENVALUES. 

REAL  A (10, 10) ,B(10, 10) ,BETA ( 10) , WK (888) , S (21 , 21 ) , M (21 , 21 ) 

REAL  G1 (3,5) , G2 (3, 5) , SUM1 , SUM2, TEMPI , TEMP2 

REAL  BZG1 (21,5) ,BZG2(21,5) , BZG1H (21 , 21 ) , BZG2H (21 , 21 ) 

REAL  BZ (21,3) ,H(5,5) ,Z1, Z2, Z3, Z4, R, L, Y, PHI (21 ) , PHI 1 (21 ) 

REAL  COMP 1 , C0MP2 

COMPLEX  ALFA(IO) , Z (10, 10) ,EIGENV(10) ,0MEGA(21) ,TEMP(21) 

COMPLEX  VECT (21, 21) .LIGVCT (10, 10) , VECTL ( 10, 10) 

INTEGER  I, J, IA, IB,N, IJOB, IZ, IER, Q, X,K,P 
INTEGER  PI ,Q2 

REAL  ZZ (21,21) ,ATRANS(10, 10) , BTRANS ( 10, 10) 

FORMAT (El 6. 8) 

FORMAT (8 (E16. 8) ) 

Q*2*P+1 

02-2*0 

NOTE  <3  IS  THE  DIMENSION  OF  THE  INPUT  MASS  AND  STIFFNESS 
MATRIX  AND  Q2  IS  THE  DIMENSION  OF  THE  OUTPUT  EIGENVECTOR 

AND  FREQUENCY  MATRIX. 

BZ (1,1) *1.0 
BZ  (1,2) *2.0 
BZ ( 1 ,3) *2.0 
Zl-L/2.0 
Z2—L/2. 0 
Z3-L/2.0 
Z4-L 
R-Zl 


CALL  CALC (R,Y,L,P, PHI, PHI 1) 
Pl-P+1 
DO  21  1*1, P 

H (2, 1+1) «PHI ( I ) 

H(4, I+PD-PHI (I) 

BZ ( 1+1 , 2) »2*PHI 1(1) 

BZ ( I+Pl , 1 ) *0. 0 
BZ ( I+Pt , 2) *0. 0 
BZ(I+1,  D-0.0 
BZ ( 1+1,3) *0.0 
H(l,  I+D-0.0 
H( i , I+Pl > -0.0 
H(2, 1+PD-0.0 


H  ( 3  ,  I +P i ) — O . 0 
H(4,I+1)«0.0 
H<5,  I+D-0.0 
CONTINUE 

H<1,1>-1.0 
H<2,1)«0.0 
H<3, l)-0.0 
H(4,l>-0.0 
H (5, 1 ) *0. 0 
R»Z2 
Y-Z4 

CALL  CALC (R, Y, L,P, PHI,  PHI 1) 

DO  30  I»1,P 

H<3,I+1)«PHI (I) 

H (3, 1+PD-PHI  (I) 

BZ ( I+Pl ,3) *2tPHIl < I > 

CONTINUE 

READ  IN  61  AND  G2 

CALL  GAIN (61, G2) 

MULTIPLY  BZ  *  G1  AND  BZ  *  G2 

DO  50  I»1,Q 
DO  60  J-1,5 
SUM 1-0. O 
SUM2-0.0 
DO  70  K*1 , 3 

TEMP1»BZ<I,K)*G1 (K, J ) 
TEMP2-BZ (I,K) *G2  <  K , J  > 

SUM1-SUM1+TEMP1 

SUM2-SUM2+TEMP2 

CONTINUE 

BZG1 <I,J)«SUM1 
BZG2(I, J)-SUM2 

CONTINUE 

CONTINUE 

MULTIPLY  BZG1*H  AND  BZG2*H 
THIS  ASSUMES  H1-H2 
DO  80  I-1,Q 
DO  90  J-1,Q 
SUM 1-0.0 
SUM2-0.0 
DO  93  K»t,5 

TEMP1-BZG1 (I,K) *H(K, J) 
TEMP2-B  ZG2  < I , K  >  *H  <  K , J ) 
SUM 1 -SUM 1 +TEMP 1 
SUM2-SUM2+TEMP2 
BZ61H ( I , J ) -SUM1 
BZG2H ( I , J ) -SUM2 

CONTINUE 


90  CONTINUE 

BO  CONTINUE 

C 

C  FORM  A  AND  B  MATRICES 
C 

DO  100  I-1,Q 

DO  110  J«1,Q 
A  < I , J) -0. 0 
Ad,  J+Q>  —  Md,J> 

Ad+Q,  J) » (S  (I , J) +BZG1H  (I ,  J ) ) 

A  < I+Q, J+Q) -BZG2H (I, J) 

B ( I , J ) — M < I , J ) 

B  < I+Q, J+Q) — — M  < I , J  > 

Bd+Q,  J) *0.0 
Bd,  J+Q) =0.0 
110  CONTINUE 

100  CONTINUE 
C 

C  SAVE  A  MATRIX  IN  FILE  AMAT  AND  B  IN  FILE  BMAT 
C  TO  PRINT  OR  SAVE  REMOVE  THE  IF  STATEMENT 

C 

IF  < I . EQ. Ill) THEN 

C  CHECK  THAT  MATRICES  WERE  READ  IN  PROPERLY 
PRINT*, '  THIS  IS  THE  A  MATRIX' 

DO  500  I»1,Q 

WRITE  (*,  6)  ( A  ( I ,  J )  ,  J**l, Q) 

500  CONTINUE 

PRINT*,'  ' 

PRINT*, 'THE  B  MATRIX  IS:’ 

PRINT*,'  ' 

DO  520  I-1,Q 

WRITE («, 6) (B(I, J) , J-1,Q> 

520  CONTINUE 
C 

ENDIF 


IA-10 
IB-10 
IZ-10 
I JOB-2 
N-Q2 
C 

CALL  EISZF  (A, IA,B, IB,N, I JOB, ALFA, BETA, Z, IZ,WK, IER> 
C 

DO  300  1-1,10 
DO  310  J-1,10 

ATRANS  < I , J ) -A ( J , I ) 

BTRANS  d,J)-B(J,I) 

310  CONTINUE 

300  CONTINUE 

CALL  EIGZF (ATRANS, I A, BTRANS, IB,N, I JOB, ALFA, BETA, 


1  LIGVCT, IZ,WK, IER) 

Q-10 

C 

DO  220  1  =1,  Q 

IF  (BETA (I) . NE.O)  THEN 
EIGENV ( I ) -ALFA ( I ) /BETA ( I > 

ELSE 

EIGENV  < I ) -999999.9 
END  IF 

220  CONTINUE 
C 

C  SORT  THE  EIGENVALUES  AND  EIGENVECTORS 
C 

X-Q-l 

221  IF  (X.NE.O)  THEN 
DO  222  I«1,X 

K-I  +  l 

COMP 1— REAL (EIGENV (K>  > 

C0MP2-REAL (EIGENV (I) ) 

IF  ( COMP 1 . GT . C0MP2 )  THEN 
TEMP ( I ) -EIGENV ( I ) 

EIGENV ( I ) -EIGENV (K) 

EIGENV (K) —TEMP ( I ) 

DO  223  J— 1 , Q 
VECTL ( J, I ) -LIGVCT ( J, I ) 

LIGVCT ( J, I ) —LIGVCT ( J , K) 

LIGVCT (J,K) -VECTL (J, I > 

VECT ( J, 1) — Z ( J, I ) 

Z ( J , I ) =Z ( J, K) 

Z(J,K)-VECT(J, I) 

223  CONTINUE 

END  IF 

222  CONTINUE 
X-X-l 
GOTO  221 
END  IF 

DO  225  1*1, Q 
DO  226  J— 1 , Q 

ZZ ( I , J) -REAL (Z ( I , J) ) 

226  CONTINUE 
225  CONTINUE 
8  FORMAT (2 (E16. 8) > 

IF  (I.EQ. Ill)  THEN 
DO  240  1-1 , Q 

PRINT*, ’EIGENVALUE  ’,  I,  ’  IS*’ 
PRINT*,’  ’ 

WRITE (*,3>EIGENV(I> 

PRINT*, ’  ’ 

OMEGA ( I ) -SORT (EIGENV ( I ) ) 

PRINT*, ’ITS  FREQUENCY  IS  IN  RAD/SEC 
PRINT*,’  ’ 

WRITE (*, 5) OMEGA ( I) 

PRINT*,’  ’ 


>c- 


w 


XT 


PRINT*, ’AND  ITS  EIGENVECTOR  IS: ’ 

PRINT*, ’  ’ 

WRITE (*,S) (Z (J, I) ,  J=1 , Q) 

PRINT*, ’  ’ 

PRINT*, ’ITS  LEFT  EIGENVECTOR  IS  :’ 

WRITE  <*, S)  (LIGVCT  < J, I) , J=1,Q) 

PRINT*,’ 

240  CONTINUE 
END  IF 
C 

END 

SUBROUTINE  CALC (R, Y, L, P, PHI , PHI 1 > 

REAL  R,Y,L,PHI (21) ,PHI1 (21) 

REAL  PI , C, D 
INTEGER  I,P 
PI»3. 14159265 
DO  311  1=1, P 
C=I*PI/L 
D=(-1)**(I+1> 

PHI ( I ) =1-C0S (C*Y) +. 5*D*C*C*Y*Y 
PHI 1(1) =C*SIN(C*R) +D*C*C*R 
311  CONTINUE 

END 

*************************************************************** 


-  C  THIS  PROGRAM  CREATES  GAIN1  (Gl)  AND  GAIN2  (G2> 

SUBROUTINE  GAIN (G, G2) 

REAL  Gl (3,5) ,G2(3,5) ,G(3,5) 

INTEGER  I,J 
C 

G( 1,2) =-1.08 

m  G ( 1 , 1 ) =8 . 60 

G( 1,3) =-3.09 
G(  1,4) —1.08 
G(l,5)=- 3.09 
6 (2, 1 ) =21 .8 
G (2, 2) =-l . 74 
t  G(2,3)=- 4.94 

G(2,4)=- 1.74 
G(2,5)»- 4.94 
G (3, 1 ) =21 . 9 
G(3,2)=- 1.74 
*  G(3,3)=- 4.94 

G(3,4)— 1.74 
G(3,5)=- 4.94 
C 

62(1, 1)=- .238 
G2(l,2)=- .016 
^  G2(l,3)»- .050 

G2(l,4)— .016 
G2(l,5)=- .050 
G2 (2, 1 ) =.643 


v 


87 


AD-A152  858  INCORPORATING  CONTROL  INTO  THE  OPTIMAL  STRUCTURAL  2/2 

DESIGN  OF  LARGE  FLEXIBL  (U)  AIR  FORCE  INST  OF  TECH 
URIGHT-PATTERSON  AFB  OH  SCHOOL  OF  ENG I.  . 

UNCLASSIFIED  T  V  HUCKENTHALER  DEC  84  AFIT/GA/AA/84D-7  F/G  22/2  NL 


G2 (2, 2) *. 233 
G2<2,3>*. 119 
G2<2,4)*-.015 
G2<2,5>-.028 
G2<3, 1) *.643 
G2<3,2)»-.015 
G2<3,3)-.028 
62 (3, 4) =.233 
G2<3,5>*.  119 
END 

**************************************************************** 

**************************************************************** 

SUBROUTINE  CALCBA<T,Hf L,P, IXPHI, IPHISQ, IPHI2, LIGVCT, 

1  RIGVCT , LBDA, PW3, PZ3) 

COMPLEX  LIGVCT (10, 10) , RIGVCT < 10, 10) , LBDA (10) 

REAL  PW3 (3) , PZ3 (3) , PIPHSQ 
REAL  P1M (21,211), P2M <21 , 21 ) , P3M  <21 , 21 > , PIS  <21 , 21 ) 

REAL  P2S<21,21) , P3S(21,21) , IXPHI <10) , IPHISQUO, 10) , IPHI2(10, 10) 

REAL  MD,  X, R, L, T,  H, MDI X, MDPH,  F, PF, PI ,  A,  B,  A2,  B2 

REAL  C, D,PXPHI, PHI , PHI 1 , PIPHI2, PHI 12, I2C,E 

REAL  PBZ (5,3) ,G1 <3,S) ,PBG1 <3, 5) , HZ (5, 5) , PBG1H<5, 5) 

REAL  PHIHF 

REAL  G2<3,3) ,PBG2(S,S) , PBG2H<5,5) 

REAL  P1A <10, 10) ,P2A<10, 10) ,P3A < 10, 10) , P1B < 10, 10) ,P2B<10, 10) 

REAL  P3B<10, 10) 

INTEGER  I,J,K,P,KI, K J 

C  THIS  PROGRAM  ASSUMES  BOTH  PAIRS  OF  APPENDAGES  ARE  KEPT  THE  SAME 
C  P1M < I , J)  IS  THE  PARTIAL  DERIVATIVE  OF  M  W.R.T  T 

C  P2M ( I , J )  IS  THE  PARTIAL  DERIVATIVE  OF  M  W.R.T.  H 

C  P3M < I , J )  IS  THE  PARTIAL  DERIVATIVE  OF  M  W.R.T.  L 

C  MD  IS  THE  MASS  DENSITY 
MD=5.22 
E*»  1  •  584E+09 
I2C-.0018 
R«1 

M2=. 136941 
C 

X-R+L 

P1M<1, 1)».6667*MD*H*<X*I3-R**3) 

P2M < 1 , 1 ) -. 66667*MD*T* < X**3-R**3) 

P3M<1, 1)=2*MD*T*H*X*X  *  4*M2*X 
C 

DO  110  I»1,P 
DO  120  J-1,P 

MDIX«MD*IXPHI <J) 

MDPH-MD* IPHISQ  < I , J ) 

P1M<1,J+1)-MDIX*H 
P1M<1,J+1+P)«MDIX*H 
P1M<1+J, 1)-MDIX*H 
PlM<i+J+P, 1)-MDIX»H 
P1M<I+1,J+1)-MDPH*H 


P1M ( 1+I+P, 1+J+P)»MDPH*H 
P1M  < 1+1 , J+l+P) *0.0 
P1M<I+1+P, J+D-0.0 
P2M<1, J+1)-MDIX*T 
P2M ( 1 , J+l+P) -MDIX*T 
P2M  <  1  +  J , 1)-MDIX*T 
P2M ( i+J+P, 1)«MDIX*T 
P2M <1+1, J+l ) »MDPH*T 
P2M  <I+t+P, J+l+P) -MDPH*T 
P2M < I+l+P, J+l ) -0. 0 
P2M  < I +1 , J+l+P) -0.0 
CONTINUE 
CONTINUE 
PI-3. 14159 

DO  130  1*1,  P 
DO  140  J«1,P 
B-J*PI 
A— I  *P  I 
A2-A*A 
B2«B*B 

C-<-l)**<I+l) 

D-(-l)**<J+l) 

C  PXPHI  IS  THE  PARTIAL  DERIVATIVE  OF  XPHI 

PXPHI— X— 2* (L/B2) *COS<B)-D*  <B2/L**3) *  < .25*X**4-.6667*X**3 

1  +.3*X*X)+.5*D*<B2/L*L)*<X**3-2*X*X+X)+2*(L/B2)+D*<B2/L**3) 

2  *.00333 

PHI— 1— COS <B)  +. 5*D*B2 
C  PHIHF  IS  PHI  EVALUATED  AT  L/2 

m  PHIHF— 1— COS  <B/2) +. 5*D* . 25*B2 

C 

HZ (1, 1)— 1.0 
HZ  < 1 , J+l ) —0.0 
HZ ( l, J+P+l ) —0.0 
HZ  (2,  0-0.0 

t  HZ (2, J+l ) —PHIHF 

HZ  <2, J+P+l) —0.0 
HZ <3, J+l) -PHI 
HZ  <3, J+l+P) -0.0 
HZ <3, l)-0.0 
HZ  <4, 1) —0.0 

*  HZ  <4, J+l ) -0.0 

HZ  (4, J+P+l ) —PHIHF 
HZ <5, 1) —0.0 
HZ (5, J+l ) —0.0 
HZ <5, J+P+l) -PHI 


PPHI1— D*B2/<L*L) 

IF  (I.NE.J)  THEN 

F-< 1/ <2*  <A-B) ) ) *L*SIN  < (A-B) )  +  <L/ <2* <A+B) ) ) 
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1  *SIN<<A+B)> 

PF-  <1/(2*  <A— B) ) ) *SIN <A-B)  +  < 1/ <2*  <A+B) ) >  *SIN  <A+B) 

ELSE 
F».5*L 
PF-.5 
END  IF 

PPHI 12— 2*  <  A2/L**3)  *B2*C*D 
C  PPHI 12  IS  THE  PARTIAL  OF  PHI’  SQ  W.R.T.  L 

PIPHSQ-1+.  1667*D*B2-D*  (B2/A2)  *CQS  (A)  +.  16667*C*A2-C*  (A2/B2) 
1  tCOS(B)  +. 0S*C*D*A2*B2  +  PF 
C  PIPHSQ  IS  THE  PARTIAL  OF  THE  INTEGRAL  OF  PHI  SQ  W.R.T.  L 
P3M <  1 ,  J+ 1 ) »MD*T*H*PXPH I +M2*PHI + I 2C*PPH I 1 
P3M  < 1 , J+P+l ) -P3M  < 1 , J+l ) 

P3M (J+l , 1 ) «P3M  < 1 , J+l ) 

P3M ( J+P+l , 1 ) «P3M ( 1 , J+l > 

P3M ( 1+1 , J+l ) «MD*T*H*PIPHSQ+I2C*PPHI 12 
P3M(I+P+1, J+P+l >-P3M <1+1, J+l) 

P3M(I+1, J+P+l) -0.0 
P3M(I+P+1, J+l) -0.0 
C  CONCLUDES  PARTIAL  OF  M 
C  THE  NEXT  PART  IS  PARTIAL  OF  S 

P I PH 1 2— 4* A2*B2* <F+C*D*L) /L**5+A2*B2*  <PF+C*D) /L**4 
P1S(I+1,J+1)«E*.2S*H*T*T*IPHI2(I, J) 

P1S<I+P+1, J+P+l) -PIS <1+1, J+l) 

P2S<I+1, J+1)-E*.0S333*T**3*IPHI2<I,J) 

P2S(I+1+P, J+1+P)-P2S(I+1, J+l) 

P3S  < 1  +  1 , J+l ) -£*. 08333*H*T**3*P1PHI2 
P3S < I+l+P, J+l+P) — P3S < 1+1 , J+l ) 

140  CONTINUE 

130  CONTINUE 

C 

K-2*P+1 

DO  150  I— 1 , K 
PlStl, I)-0.0 
P1S<I, l)-0.0 
P2S (1,1) —0.0 
P2S  < 1 , 1 ) -0. 0 
P3S  <  1,  I )  — 0. 0 
P3S (1,1) —0 . 0 
150  CONTINUE 

C 

DO  160  I-1,P 
DO  170  J-i,P 

P1S<I+1, J+l+P)-0.0 
PlS<I+P+l,J+l)-0.0 
P2S <1+1, J+P+l) -0.0 
P2S  < I+P+i , J+l ) -0. 0 
P3S <1+1, J+P+l) -0.0 
P3S ( I+P+l ,J+t)-0.0 
170  CONTINUE 

160  CONTINUE 

C 

DO  180  I«1,P 
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n  n 


A-I*PI 

PBZ ( I  +  i , 2) — (A/ (L*L) ) *SIN ( . 5*A> 
PBZ < I+P+l , 3)  —  <A/ <L*L>  >  *SIN (A/2. 03 
PBZ(I+1,1)«0.0 
PBZ  ( I+l+P,  1 )  “0.0 
PBZ (1+1,3) -0.0 
PBZ ( I+l+P, 2) -0.0 
180  CONTINUE 

PBZ (1, 1) -0.0 
PBZ (1,2) -0.0 
PBZ (1,3) —0.0 

MULTIPLY  PBZ  *G1 
K-2*P+1 

CALL  GAIN(G1,G2) 

CALL  MATMUL(K,3,3,PBZ,G1,PBG1) 

C  MULTIPLY  PBGltH 

CALL  MATMUL (K, 5, K, PBG1 , HZ, PBG1H) 
CALL  MATMUL (K, 3,9, PBZ, G2,PBG2) 

CALL  MATMUL (K, 5, K, PBG2, HZ, PB62H) 

C  ASSEMBLE  PARTIAL  OF  A  AND  PARTIAL  OF  B 
K»2*P+1 
DO  190  I— 1 , K 
DO  200  J— 1,K 
KJ-J+K 
KI-I+K 
P1A ( I , J) — 0.0 
P2A(I, J)» 0.0 
P3A ( I, J) — 0.0 
P1A(I,KJ)  —  (P1M(I, J)*2) 

P2A ( I , K J ) — — 2*P2M ( I , J ) 

P3A ( I , KJ ) — -2*P3M ( I , J ) 
P1A(KI,J)-2*P1S(I,J) 

P2A (K I , J ) -2*P2S ( 1 , J ) 

P3A (KI , J) - (2*P3S ( I , J) +PBG1H ( I , 
P1A (KI , KJ) — 0.0 
P2A(KI,KJ)— 0,0 
P3A (KI , KJ ) -PBG2H (I, J) 

200  CONTINUE 

190  CONTINUE 

C 

DO  210  I— 1 ,  K 
DO  220  J— l,K 
KI-I+K 
KJ-J+K 

P1B(I, J)--2*P1M(I, J) 

P2B ( I , J ) »-2*P2M ( I , J ) 

P3B(I,J)— 2*P3M(I,J) 

P1B(I,KJ)— 0.0 
P2B ( I , KJ ) —0.0 
P3B (I,KJ)-0.0 

P1B(KI, J)-0. 0 
P2B ( K I , j ) -0 . 0 


■fuV.,.v  J 


.'.••/.'A  •.  •  • 
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P3B(KIf J)-0.0 
P1B(KI,KJ)— *2*P1M(I,  J) 

P2B (KI , KJ) — 2*P2M (I, J) 

P3B (KI,KJ) — 2*P3M (I, J) 

220  CONTINUE 
210  CONTINUE 

CALL  GR ADNT <  P 1 A  , P2 A , P3 A , P 1 B , P2B , P3B , L I G  VCT , 

1  RIGVCT  ,LBDA,PH3, PZ3) 

END 

**************************************************************** 
«*************************f ************************************* 

SUBROUTINE  GRADNT(P1A,P2A,P3A,P1B,P2B,P3B,LIGVCT, 

1  RIGVCT, LBDA,PW3,PZ3) 

C  THIS  PROGRAM  CALCULATES  THE  GRADIENTS  OF  OMEGA 
C  AND  ZETA  W.R.T  THE  PARAMETERS  T,  H  ,  L 
C 

REAL  PI A (10, 10) , P2A ( 10, 10)  , P3A <10,10) , P1B ( 10, 10),P2B(10, 10) 
REAL  P3B ( 10, 10) 

COMPLEX  LIGVCT ( 10, 10) , RIGVCT (10, 10) ,LBDA(10) ,LTRANS(10, 10) 
COMPLEX  LPA1 (10,10) ,LPA1R< 10, 10) ,LPB1 (10, 10) ,LPB1R(10, 10) 
COMPLEX  LPA2 ( 10, 10) , LPA2R ( 10, 10) , LPB2 ( 10, 10) , MAT ( 10, 10) 
COMPLEX  LPB2R ( 10, 10) ,LPA3(10, 10) , LPA3R ( 10, 10) ,LPB3(10, 10) 
COMPLEX  LPB3R  < 10, 10) , PI LB DA ( 10) , P2LBDA ( 10) , P3LBDA ( 10) 
INTEGER  I,J,K,P,DIM 
REAL  PW3(3),PZ3(3),DEN0M,NUM1,NUM2 
C 

P-2 

DIM-2* <2*P+1) 

C 

C  TRANSPOSE  LIGVCT 
DO  20  I— 1 , DIM 
DO  30  J— 1,DIM 

LTRANS (I, J) -LIGVCT (J, I) 

30  CONTINUE 

20  CONTINUE 

C 

C  CALCULATE  THE  PARTIAL  DERIVATIVES  OF  LAMBDA 
C 

DIM-10 

C 

c 

CALL  MATCR (LTRANS, PI A, LPA1 ) 

CALL  MTCPLX(LPA1, RIGVCT, LPA1R) 

CALL  MATCR < LTRANS, P1B,LPB1) 

CALL  MTCPLX(LPB1, RIGVCT, LPBIR) 

C 

C  PRINT*, 'D  LAMDA  /DT  IS  ' 

DO  40  I-1,DIM 
DO  50  J-l , DIM 

MAT  < I , J) — LPA1R  < I , J ) -LBDA ( I ) *LPB1R  < I , J ) 


SO  CONTINUE 

P1LBDA(I)»MAT(I, I> 

C  PRINT*, PILBDA(I) 

40  CONTINUE 

C 

C  CALCULATE  PARTI ALS  OF  LAMBDA  W.R.T.  P2 
CALL  MATCR ( LTR ANS , P2A , LPA2 > 

CALL  MTCPLX (LPA2, RIGVCT, LPA2R) 

CALL  MATCR ( LTR ANS , P2B , LPB2 > 

CALL  MTCPLX (LPB2, RIGVCT, LPB2R) 

C 

C  PRINT*, ’D  LAMDA/DH  IS  ' 

DO  60  1*1 , DIM 
DO  70  J*1,DIM 

MAT ( I , J ) *LPA2R ( I , J ) — LBDA ( I ) «LPB2R ( I , J ) 

70  CONTINUE 

P2LBDA ( I ) »MAT ( I ,  I> 

C  PRINT* , P2LBDA  < I ) 

60  CONTINUE 

C 

C  CALCULATE  PART I ALS  OF  LAMBDA  W.R.T.  P3,  P3-L 
C 

CALL  MATCR (LTRANS, P3A, LPA3) 

.  CALL  MTCPLX (LPA3, RIGVCT, LPA3R) 

CALL  MATCR (LTRANS, P3B,LPB3) 

CALL  MTCPLX (LPB3, RIGVCT, LPB3R) 

C 

C  PRINT*, »D  LAMBDA  /DL  IS  ' 

DO  80  I-t,DIM 
DO  90  J-l.DIM 

MAT (I, J) -LPA3R (I, J) -LBDA  < I ) *LPB3R  < I , J ) 

90  CONTINUE 

P3LBDA < I ) *»MAT ( I ,  I) 

C  PRINT* , P3LBDA  < I ) 

80  CONTINUE 

C 

C  FORM  GRADIENTS  OF  W1  AND  ZETA1 
C  PW3 ( I )  IS  THE  PARTIAL  OF  W(I)  W.R.T.  PI 
C  PZ3 (I)  IS  THE  PARTIAL  OF  ZETA(l)  W.R.T.  PI 
C 

c 

DENOM*SQRT ( (REAL (LBDA (3) ) ) **2+ (AIMAG (LBDA (3) ) ) **2> 

C 

PW3 ( 1 ) * (REAL (P1LBDA (3) ) +AIMAG (P1LBDA (3) ) ) /DENOM 
PW3 (2) »(REAL (P2LBDA (3) ) +AIMAG (P2LBDA (3) ) ) /DENOM 
PW3(3)»(REAL(P3LBDA(3) > +AIMAG (P3LBDA(3> ) ) /DENOM 
C 

NUM1- (REAL (LBDA <3> ) )**2 

NUM2- (REAL (LBDA (3) ) ) * (AIMAG (LBDA (3) ) ) 

PZ3(t>—  REAL (P1LBDA (3) )  /DENOM  -MNUM1*  (REAL (P1LBDA (3)  )  )+NUM2 
1  «AIMAG (P1LBDA (3) ) ) /DENOM* *3 


PZ3<2)»-REAL<P2LBDA<3>  > /DENOM  ♦  (NUhl*  (REAL <P2LBDA<3> )  )  +NUM2 
1  tAIMAG (P2LBDA  <3) ) ) /DEN0M**3 
C 

PZ3  <3) *<— REAL  (P3LBDA  (3) )  /DENOh)  ♦  (NUM1* <REAL<P3LBDA< 3> ) )♦ 

1  NUM2 tAIMAG  (P3LBDA  (3) ) )  /DENOMt*3 
C 

END 

SUBROUTINE  MATMUL < NR t , NR2, NCOL2, A, B, C> 

INTEGER  NR1,NR2, NC0L2, I,  J,K 
REAL  SUM1, TEMPI 

REAL  A(NR1,NR2) , B <NR2, NC0L2) , C  <NR1 , NC0L2) 

DO  10  I-1,NR1 
DO  20  J»1 , NC0L2 
SUM 1-0.0 

DO  30  K— 1 , NR2 
TEMP1-A<I,K) *B<K, J) 

SUM 1 —SUM 1 ♦TEMPI 
C<I, J)-SUM1 
30  CONTINUE 

20  CONTINUE 

10  CONTINUE 

END 
C 

SUBROUTINE  MTCPLX (A, B, C) 

C  THIS  SUBROUTINE  MULTIPLIES  TWO  COMPLEX  MATRICES  TOGETHER 
C 


COMPLEX  C<10,10) 

INTEGER  NR1 y  NR2, NC0L2 

COMPLEX  A <10, 10) ,B(10y 10) 
COMPLEX  SUMt, TEMPI 
INTEGER  I , J, K 
C 

NR1-10 

NR2-10 

NC0L2-10 

C  MULTIPLY  GIVEN  MATRICES 
C 

DO  10  1*1, NR1 
DO  20  J*1 , NC0L2 
SUMisao 

DO  30  K-1,NR2 

TEMP1*A(I,K) *B  <K, J> 

SUM 1*SUM1 ♦TEMPI 
30  CONTINUE 

C<I,J)*SUM1 
20  CONTINUE 

10  CONTINUE 
END 
C 
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C  THIS  SUBROUTINE  MULTIPLIES  A  COMPLEX  AND  A  REAL  MATRIX 
C  TOGETHER  ,  THE  MATRICES  ARE  10  X  10 
C 

SUBROUTINE  MATCR<A,B,C) 

INTEGER  NR 1 , NR2 , NC0L2 
COMPLEX  A<10,10>,C<10,10> 

REAL  B ( 10, 10) 

COMPLEX  SUM1, TEMPI 
INTEGER  I,J,K 
C 

NR1-10 

NR2-10 

NC0L2-10 

C  MULTIPLY  GIVEN  MATRICES 
C 

DO  10  I-1,NR1 
DO  20  J«1,NC0L2 
SUM1-0.0 
DO  30  K«1,NR2 

TEMP1*A<I,K)*B<K,J> 

SUM1-SUM1+TEMP1 
30  CONTINUE 

C(I, J)-SUM1 
20  CONTINUE 

10  CONTINUE 

END 

*******t**t***t**t*tt**tt**t*tt*t**«*tttttt*t****s***t***t*t*** 

*******t****tt****ttt*ttt**tt**U**tt*tt*t*t**ttt ************** 


CONMIN  OUTPUT  FOR  CASE  2 

This  appendix  contains  a  more  detailed  listing  of  the 
output  referred  to  in  Chapter  4,  Case  2.  The  listings 
provide  the  initial  values  for  the  decision  variables,  t,  h, 
and  L  and  the  performance  criteria  used  for  the  particular 
run,  i.e.  DELFUN,  DAB  FUN,  CT,  and  THETA.  The  final 
optimization  information  is  also  shown  with  the  reason  for 
termination  and  the  number  of  iterations  required  for 
convergence  or  the  number  of  iterations  performed  without 


convergence. 


C  Q  N  M  I  N 
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FORTRAN  PROGRAM  FOR 
CONSTRAINED  FUNCTION  MINIMIZATION 
********************* 


* 

* 

* 

* 

* 

* 


CONSTRAINED  FUNCTION  MINIMIZATION 


CONTROL  PARAMETERS 


I PR I NT  NDV 

ITMAX  NCON 

NS IDE  ICNDIR 

NSCAL  NFDG 

2  3 

40 

2 

6  4 

0  1 

LINOBJ  ITRM 

N1 

N2 

N3  N4 

NS 

0  3 

3 

8 

51  31 

102 

CT 

CTMIN 

CTL 

CTLMIN 

10000E+00 

• 

40000E— 02 

-. 10000E-01 

. 1 OOOOE— 02 

THETA 

PHI 

DELFUN 

DABFUN 

. 1 OOOOE+O 1 

• 

30000E-K>1 

. 1 OOOOE— 03 

- 40507E-03 

FDCH 

FDCHM 

ALPHAX 

AB0BJ1 

. lOOOOE— 01 

• 

lOOOOE— 01 

. 10000E+00 

. 10000E+00 

LOWER  BOUNDS  ON  DECISION  VARIABLES  (VLB) 

l)  . lOOOOE— 03  . 10000E— 05  . 10000E-03 

UPPER  BOUNDS  ON  DECISION  VARIABLES  (VUB) 

1)  .  10000E+U  . 10000E+1 1  . lOOOOE+ll 

ALL  CONSTRAINTS  ARE  NON-LINEAR 


INITIAL  FUNCTION  INFORMATION 

OBJ  »  . 405072E+0O 

DECISION  VARIABLES  (X-VECTOR; 

1)  . 97000E-02  . 50000E+00  .40000E+01 

CONSTRAINT  VALUES  (G-VECTOR; 

1 )  . 42565E+00  - . 33929E-02 


ITER  »  i  OBJ  -  .40507E+00  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 97000E— 02  . 50000E+00  .40000E+01 

CONSTRAINT  VALUES  (G- VECTOR; 

. 42365E+00  - - 35929E-02 


ITER  -  2  OBJ  -  . 40507E+00  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 97000E-02  . 30000E+00  .40000E+01 

CONSTRAINT  VALUES  <G- VECTOR/ 

1 )  . 42363E+00  - . 35929E-02 


ITER  -  3  OBJ  *  .40507E+00  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  (X-VECTOR) 

1)  .  97000E— 02  .30000E+00  . 40000E+01 

CONSTRAINT  VALUES  <G- VECTOR) 

1  >  .  42365E-K>0  - .  35929E-02 


ITER  -  4  OBJ  *  . 40507E+00  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  (X-VECTOR) 

i)  . 97000E— 02  . 50000E+00  .40000E+01 

CONSTRAINT  VALUES  (G- VECTOR) 

1 )  . 42565E+00  -. 35929E-02 


ITER  =  3  OBJ  -  . 41520E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 99422E-02  .30001E+00  .40001E+01 

CONSTRAINT  VALUES  (G- VECTOR) 

1 )  . 28323E+00  - . 23034E-02 


ITER  -  b  OBJ  -  . 43479E+00 


DECISION  VARIABLES  (X-VECTOR) 

1)  . 1041 IE-01  . 30002E+00  . 40002E+01 

CONSTRAINT  VALUES  (G-VECTOR) 


1 )  . 3987 1 E-02  - . 28098E-04 

ITER  »  7  OBJ  *  . 43479E+00  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  (X-VECTOR/ 

1)  .1041 IE-01  . 50002E+00  .40002E+01 

CONSTRAINT  VALUES  (G- VECTOR) 

1 )  . 3987 1 E-02  - . 28098E-04 

ITER  «  8  OBJ  »  . 43479E+00  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 1041 IE-01  . 50002E+00  .40002E+01 

CONSTRAINT  VALUES  <G- VECTOR) 

1 )  . 3987 l E-02  - . 28098E— 04 

FINAL  OPTIMIZATION  INFORMATION 

OBJ  *  . 434785E+0G 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 1041 IE— 01  .50002E+00  .40002E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  .3987 IE— 02  — . 28098E-04 

THERE  ARE  2  ACTIVE  CONSTRAINTS 
CONSTRAINT  NUMBERS  ARE 
1  2 

THERE  ARE  O  VIOLATED  CONSTRAINTS 
THERE  ARE  *  0  ACTIVE  SIDE  CONSTRAINTS 
TERMINATION  CRITERION 

ABS ( 1-OBJ ( I— 1 ) /OBJ ( I ) )  LESS  THAN  DELFUN  FOR  3  ITERATIONS 
NUMBER  OF  ITERATIONS  -  8 

OBJECTIVE  FUNCTION  WAS  EVALUATED  24  TIMES 

CONSTRAINT  FUNCTIONS  WERE  EVALUATED  24  TIMES 

GRADIENT  OF  OBJECTIVE  WAS  CALCULATED  5  TIMES 

GRADIENTS  OF  CONSTRAINTS  WERE  CALCULATED  S  TIMES 


*  * 

*  FORTRAN  PROGRAM  FOR  * 

*  * 

*  CONSTRAINED  FUNCTION  MINIMIZATION  * 

*  * 


ft************************** 

CONSTRAINED  FUNCTION  MINIMIZATION 


CONTROL  PARAMETERS 


I PRINT 

2 

NDV 

3 

ITMAX  NCON 

40  2 

NS IDE  ICNDIR 

6  4 

NSCAL  NFDG 

0  1 

LINOBJ 

0 

ITRM 

3 

N1  N2 

S  8 

N3  N4 

51  51 

N5 

102 

CT 

10000E+01 

CTMIN 
. 40000E-02 

CTL 

-. lOOOOE— 01 

CTLMIN 
. 1 OOOOE— 02 

THETA 
. 20000E+0 1 

PHI 

.50000E+01 

DELFUN 
. 1 OOOOE— 03 

DABFUN 
. lOOOOE— 01 

FDCH 

. lOOOOE— 01 

FDCHM 
. lOOOOE— 01 

ALPHA X 
. 1 OOOOE+OO 

AB0BJ1 
. 1 OOOOE+OO 

LOWER  BOUNDS  ON  DECISION  VARIABLES  (VLB) 

t>  . lOOOOE— 05  • lOOOOE— 05  . 10000E-05 

UPPER  BOUNDS  ON  DECISION  VARIABLES  (VUB) 

1)  . 10000E+1 1  . 10000E+1 1  . 1 0000E+ 1 1 

ALL  CONSTRAINTS  ARE  NON-LINEAR 


INITIAL  FUNCTION  INFORMATION 

OBJ  -  .501 120E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  .  I5000E-01  •  40000E+00  .40000E+01 

CONSTRAINT  VALUES  (G-VECTOR) 
i )  - . 23225E+0 1  . i 0955E-0 1 


ITER  »  1  OBJ  -  .501 12E+00  NO  CHANGE  IN  OBJ 


40000E+01 


DECISION  VARIABLES  (X-VECTOR) 

1 >  .  1 5000E-0 1  . 40000E+00 

CONSTRAINT  VALUES  (G-VECTOR) 

23223E+01  . 10953E-01 

ITER  *  2  OBJ  *  .5011 2E+00  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  (X- VECTOR) 

1)  . 15000E-01  . 40000E+00  . 40000E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1 )  - . 23223E+0 1  . 1 0955E-0 1 


ITER  -  3  OBJ  -  . 49979E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 14959E— 01  . 40000E+00  . 40002E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  -.22962E+01  . 10S71E— 01 


ITER  *  4  OBJ  -  . 49846E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  • 14919E— Ol  . 40000E+00  . 40004E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  -.22700E+01  .  10787E— 01 


ITER  »  5  OBJ  -  •  497 1 4E+00 

DECISION  VARIABLES  (X-VECTOR) 

l)  » 14879E-01  . 40000E+00  .4C006E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  — . 22440E+01  . 10702E— 01 


ITER  ■  6  ■  OBJ  -  . 49582E+00 

DECISION  VARIABLES  (X-VECTOR) 

. 14938E-01  . 40000E+00  . 40008E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  -.221S0E+01  . 10617E-01 


ITER 


7 


OBJ 


49450E+00 


DECISION  VARIABLES  (X- VECTOR) 

1)  . 14798E-01  • 40000E+00 


. 40010E+01 


CONSTRAINT  VALUES  <G- VECTOR) 

i)  -.21921E+01  . 10532E-01 


ITER 


.49319E+00 


DECISION  VARIABLES  (X-VECTOR) 

l)  ► 14758E-01  . 40000E+00  .40012E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  -.21663E+01  . 10446E— 01 


ITER 


OBJ  - 


. 49 1 88E+00 


DECISION  VARIABLES  (X-VECTOR) 

1)  . 14718E— 01  . 40000E+00  .40014E+01 

CONSTRAINT  VALUES  <G- VECTOR) 

1)  -.2I406E+01  . 10360E— 01 


ITER 


. 49188E+00 


NO  CHANGE  IN  OBJ 


DECISION  VARIABLES  (X-VECTOR) 

1)  • 147I8E— 01  . 40000E+00  .40014E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  -.2140AE+01  . 10360E-01 


FINAL  OPTIMIZATION  INFORMATION 


.491878E+00 


DECISION  VARIABLES  (X-VECTOR) 

1)  . 14718E-01  . 40000E+00  .40014E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

l)  -. 21406E+01  .103A0E-01 

THERE  ARE  O  ACTIVE  CONSTRAINTS 

THERE  ARE  1  VIOLATED  CONSTRAINTS 

CONSTRAINT  NUMBERS  ARE 
2 


THERE  ARE 


O  ACTIVE  SIDE  CONSTRAINTS 


■  '*'■  . 


-•  v  vrsBSfigft 


TERMINATION  CRITERION 

TEN  CONSECUTIVE  ITERATIONS  FAILED  TO  PRODUCE  A  FEASIBLE  DESIGN 
NUMBER  OF  ITERATIONS  -  10 

OBJECTIVE  FUNCTION  WAS  EVALUATED  31  TIMES 

CONSTRAINT  FUNCTIONS  WERE  EVALUATED  31  TIMES 

GRADIENT  OF  OBJECTIVE  WAS  CALCULATED  9  TIMES 

GRADIENTS  OF  CONSTRAINTS  WERE  CALCULATED  9  TIMES 


CONMIN  ODTPDT  fOR  £ASE  £ 

Case  3  elliaimted  the  equality  constraints  of  Case  2 
and  replaced  then  with  inequality  constraints.  A  region  was 
specified  to  aaintain  the  natural  frequency  of  the  third 
eigenvalue,  in  order  to  provide  CONMIN  sore  flexibility  in 
the  optimization  process.  Only  a  few  of  the  many  trials  are 
listed  to  demonstrate  that  more  than  one  local  minimum  exists 
due  to  the  linearity  of  the  objective  function  and  the 
gradients  of  the  eigenvalues  with  respect  to  the  structural 


parameters . 


*  CONMIN  * 

*  * 

*  FORTRAN  PROGRAM  FOR  * 

•  *  * 

*  CONSTRAINED  FUNCTION  MINIMIZATION  * 

*  * 
*************************** 

CONSTRAINED  FUNCTION  MINIMIZATION 

* 

CONTROL  PARAMETERS 


IPRINT 

NDV 

ITMAX 

NCON 

NSIDE 

ICNDIR 

NSCAL 

2 

3 

40 

2 

6 

4 

0 

L I NOB J 

ITRM 

N1 

N2 

N3 

N4 

N5 

0 

3 

S 

8 

51 

51 

102 

20000E+01 


CTMIN  CTL 

. 40000E— 02  10000E-01 


CTLMIN 


. 10000E-02 


THETA 
. 30000E+01 

FDCH 

. 10000E— 01 


PHI 

•50000E+01 

FDCHM 
. 10000E-01 


DELFUN 
. 10000E-03 

ALPHA X 
. 10000E+00 


LOWER  BOUNDS  ON  DECISION  VARIABLES  (VLB) 

1)  . 10000E— 03  . 10000E-05  . 10000E-05 

UPPER  BOUNDS  ON  DECISION  VARIABLES  <VUB> 

1 )  . 10000E+1 1  . 10000E+1 1  . 10000E+1 l 

ALL  CONSTRAINTS  ARE  NON-LINEAR 


DABFUN 

10000E+00 

AB0BJ1 

10000E+00 


INITIAL  FUNCTION  INFORMATION 

OBJ  •  . 402372E+00 

DECISION  VARIABLES  (X-VECTOR) 
l)  . 10417E-01  . 50000E+00 

CONSTRAINT  VALUES  (G- VECTOR) 

1)  .36013E+00  76013E+00 


•37000E+01 


ITER  -  1  OBJ  *  .40237E+00  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 10417E— 01  .30000E+00  .37000E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

36013E+00  -.76013E+00 


ITER  *  2  OBJ  -  .40237E+00  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  (X -VECTOR) 

i)  . 10417E-01  . 50000E+00  .37000E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1 )  . 360 1 3E+00  -.7601 3E+00 


ITER  *  3  OBJ  -  . 39224E+00 

DECISION  VARIABLES  < X-VECTOR) 

1)  . 10491E— 01  .48477E+00  .36938E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  .33718E+00  -. 75718E+00 


ITER  -  4  OBJ  »  . 19333E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 1 1930E— 01  . 21629E+00  .33882E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  — . 26099E+00  -. 13901E+00 


ITER  -  3  OBJ  -  . 19333E+00  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  (X-VECTOR) 

l)  . 1 1930E-01  . 21629E+00  .33882E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

l )  - . 26099E+00  - . 1 390 l E+00 


ITER  -  6  OBJ  -  . 18221 E+00 


DECISION  VARIABLES  (X-VECTOR) 

1)  . 12963E-01  • 18786E+00  .35831E+01 


CONSTRAINT  VALUES  (G-VECTOR) 


1>  —.9171 IE-04  -.39991E+00 

ITER  -  7  OBJ  -  .171 26E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 121S6E-01  .18784E+00  .35831E+01 

CONSTRAINT  VALUES  (G- VECTOR) 

1 )  - . 40000E+00  - . 2863 1 E-06 

FINAL  OPTIMIZATION  INFORMATION 

OBJ  -  . 171257E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 12186E-01  . 18784E+00  .33831E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  — . 40000E+00  — . 28651 E-06 

THERE  ARE  1  ACTIVE  CONSTRAINTS 
CONSTRAINT  NUMBERS  ARE 
2 

THERE  ARE  O  VIOLATED  CONSTRAINTS 

THERE  ARE  0  ACTIVE  SIDE  CONSTRAINTS 

TERMINATION  CRITERION 

ABS (OBJ ( I ) —OBJ ( I— 1 ) )  LESS  THAN  DABFUN  FOR  3 


NUMBER  OF  ITERATIONS  - 

7 

OBJECTIVE  FUNCTION  WAS  EVALUATED 

19 

TIMES 

CONSTRAINT  FUNCTIONS  WERE 

EVALUATED 

19 

TIMES 

GRADIENT  OF  OBJECTIVE  WAS 

CALCULATED 

5 

TIMES 

ITERATIONS 


GRADIENTS  OF  CONSTRAINTS  W'  RE  CALCULATED 


3  TIMES 


************************* 

* 

*  C  0  N  M  I  N 

* 

*  FORTRAN  PROGRAM  FOR 

* 

*  CONSTRAINED  FUNCTION  MINIMIZATION 

* 

41************************ 


*  4t 
* 
* 
* 
* 
* 
* 
* 

*  * 


CONSTRAINED  FUNCTION  MINIMIZATION 
CONTROL  PARAMETERS 


IPRINT 

2 

NDV 

3 

ITMAX  NCON 

40  2 

NS  IDE  ICNDIR 

6  4 

NSCAL  NFDG 

0  l 

LINOBJ 

0 

ITRM 

3 

N1  N2 

5  a 

N3  N4 

31  51 

N5 

102 

CT 

1C000E+00 

CTMIN 
. 40000E-02 

CTL 

10000E— 01 

CTLMIN 
. 10000E-02 

THETA 
. 10000E+01 

PHI 

. 50000E+0 1 

DELFUN 

.10000E-03 

DABFUN 
. 43430E— 03 

FDCH 

. 10000E— 01 

FDCHM 
. lOOOOE-Ol 

ALPHA X 
. 10000E+00 

AB0BJ1 
. 10000E+00 

LOWER  BOUNDS  ON  DECISION  VARIABLES  (VLB) 

1)  . 10000E-05  . 10000E-03  . 10000E-05 


UPPER  BOUNDS  ON  DECISION  VARIABLES  <VUB> 

l)  .  IOOOOE+-1 1  .  iOOOOE+t 1  . IOOOOE+11 


ALL  CONSTRAINTS  ARE  NON-LINEAR 


INITIAL  FUNCTION  INFORMATION 

OBJ  »  . 434304E+00 

DECISION  VARIABLES  (X-VECTOR) 

i>  . 10400E-01  . 50000E+00  . 40000E+01 

CONSTRAINT  VALUES  (G- VECTOR) 

i)  — . 22037E+00  17963E+00 


ITER  -  l  OBJ  ■  .42158E+00 
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40000E+0 1 


DECISION  VARIABLES  (X-VECTOR) 

1)  . 10093E— 01  .49999E+00 

CONSTRAINT  VALUES  (G-VECTOR) 

i )  -. 40260E+00  . 2600 IE-02 


ITER  «  2  OBJ  -  . 38001 E+OO 

DECISION  VARIABLES  (X-VECTOR) 

1 )  .1 1363E-0 1  . 40405E+00 

CONSTRAINT  VALUES  <G- VECTOR) 

1 )  - . 73960E— 03  - . 39926E+00 


ITER  -  3  OBJ  *  • 35673E+00 

DECISION  VARIABLES  (X-VECTOR) 
i )  . 10668E-0 1  . 40403E+00 

CONSTRAINT  VALUES  (G- VECTOR) 

1)  — . 40023E+00  . 22547E— 03 


ITER  »  4  OBJ  -  .31854E+00 

DECISION  VARIABLES  (X-VECTOR) 

1 )  . 12125E-0 1  . 3 1952E+00 

CONSTRAINT  VALUES  (G-VECTOR) 

1 )  67936E— 03  - . 39932E+00 


ITER  -  5  OBJ  -  .31851E+00 

DECISION  VARIABLES  (X-VECTOR) 
i)  . 12127E-01  . 31943E+00 

CONSTRAINT  VALUES  (G-VECTOR) 

1 )  - . 1 3330E-04  - . 39999E+00 


ITER  »  6  OBJ  -  . 31851E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 12127E-01  . 31943E+00 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  —.6221 IE— 07  -.40000E+00 


39638E+01 


39638E+01 


39380E+01 


39380E+01 


39380E+0 1 


ITER 


7 


OBJ 


31833E+00 


39379E+01 


DECISION  VARIABLES  (X-VECTOR) 
i)  . 12138E-01  .31 899E+00 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  . 32826E-02  40328E+00 

FINAL  OPTIMIZATION  INFORMATION 

OBJ  *  .318350E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 12138E-01  . 31899E+00  .39379E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1 )  . 32826E— 02  - . 40328E+00 

THERE  ARE  1  ACTIVE  CONSTRAINTS 
CONSTRAINT  NUMBERS  ARE 
1 

THERE  ARE  O  VIOLATED  CONSTRAINTS 

THERE  ARE  O  ACTIVE  SIDE  CONSTRAINTS 

TERMINATION  CRITERION 

ABS ( OBJ ( I ) —OBJ ( I — 1 ) )  LESS  THAN  DABFUN  FOR  3  ITERATIONS 

NUMBER  OF  ITERATIONS  «  7 

OBJECTIVE  FUNCTION  WAS  EVALUATED  22  TIMES 

CONSTRAINT  FUNCTIONS  WERE  EVALUATED  22  TIMES 

GRADIENT  OF  OBJECTIVE  WAS  CALCULATED  7  TIMES 

GRADIENTS  OF  CONSTRAINTS  WERE  CALCULATED  7  TIMES 


APPENDIX  F 


CONMIN  OUTPUT  FOR  CASE  £ 

In  Case  4  the  frequency  for  the  third  eigenvalue  was 
constrained  to  remain  between  4.2  and  5.2  and  the  first 
eigenvalue  was  constrained  to  the  left  half  plane.  In  Case  3 
the  objective  funciton  was  minimized  and  the  frequency  of  the 
third  eigenvalue  maintained  in  a  specified  region,  but  the 

l 

first  eigenvalue  was  unconstrained  and  hence  moved  to  an 
unstable  region.  Since  this  instability  can  occur,  the  first 
eigenvalue  was  constrained  to  the  left  half  plane.  Three  of 
the  many  trials  attempted  are  listed  here  as  referenced  in 
Case  4  . 


Ill 


* 

* 

* 

* 

* 

* 


CONSTRAINED  FUNCTION  MINIMIZATION 


* 

* 

* 

* 

* 

* 


*************************** 


CONSTRAINED  FUNCTION  MINIMIZATION 


CONTROL  PARAMETERS 


I PRINT  NDV 

ITMAX  NCON 

NS IDE  ICNDIR 

NSCAL  NFDG 

2  3 

40 

3 

6  '  4 

0  1 

LINOBJ  ITRM 

N1 

N2 

N3  N4 

N5 

0  3 

3 

S 

51  51 

102 

CT 

CTMIN 

CTL 

CTLMIN 

10000E+01 

■ 

40000E-02 

10000E-01 

. 10000E-02 

THETA 

PHI 

DELFUN 

DABFUN 

. 30000E+0 1 

• 

50000E+01 

. 10000E-03 

. 43500E-03 

FDCH 

FDCHM 

ALPHA X 

ABOB J 1 

. lOOOOE— 01 

• 

10000E-01 

. 10000E+00 

. 10000E+00 

LOWER  BOUNDS  ON  DECISION  VARIABLES  (VLB) 

l)  . lOOOOE— 05  . 10000E-05  . 10000E-05 


UPPER  BOUNDS  ON  DECISION  VARIABLES  (VUB) 

l)  . 10000E+1 1  . 10000E+1 1  . 10000E+1 1 

ALL  CONSTRAINTS  ARE  NON-LINEAR 


INITIAL  FUNCTION  INFORMATION 

OBJ  =  . 434997E+00 

DECISION  VARIABLES  (X-VECTOR) 

i)  . 10417E-01  .50000E+00  .40000E+01 

CONSTRAINT  VALUES  (G- VECTOR) 

1)  — .  81038E+00  18962E+00  -.S0229E-01 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  -  1  OBJ  *  . 43500E+00  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  <X-VECTOR) 

1)  . 10417E-01  .50000E+00  .40000E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  -. 8103SE+00  -. 18962E+00  -.80229E-01 


ITER  *  2  OBJ  *  . 43340E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 10378E-01  . 50000E+00  . 40000E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

i>  -.83335E+00  -. 16665E+00  -.80799E-01 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  =  3  OBJ  *  . 42943E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  .  10283E— 01  . 50000E+00  .40000E-K>1 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  — . 89042E+00  -. 10958E+00  -.77555E-01 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  *  4  OBJ  »  . 42755E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 10238E— 01  . 50000E+00  . 40000E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

i)  — . 91730E+00  — . 82699E-01  -. 33336E-01 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  *  5  OBJ  -  . 42737E+00 

DECISION  VARIABLES  (X-VECTOR) 

i)  . 10234E-01  . 50000E+00  .40000E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

i)  — . 91989E+00  -. 80109E-01  -.44576E-02 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  - 


6 


OBJ  * 


42037E+00 


39936E+01 


DECISION  VARIABLES  (X-VECTOR) 
i)  . 10440E— 01  . 482B7E+00 

CONSTRAINT  VALUES  (G-VECTOR) 

i)  -.S3074E+00  -.14926E+00  -.81351E-01 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  -  7  OBJ  *  .41684E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 103S2E-01  . 48287E+00  .39936E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  — . 90293E+00  -.97074E— 01  -.21008E-01 

*#  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  *  8  OBJ  *  . 41684E+00  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 10352E— 01  . 48287E+00  .39936E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

i)  -.90293E+00  — . 97074E— 01  -.21008E-01 


ITER  »  9  OBJ  »  .41674E+00 

DECISION  VARIABLES  (X-VECTOR) 

l)  . 10350E— 01  . 48287E+00  .39936E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  — • 90438E+00  - . 93625E-0 1  -. 20532E-02 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  *  10  OBJ  »  .41 673E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 10350E-01  . 482B7E+00  .39936E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

X)  -. 90452E+00  -.95484E-01  -.67406E-05 


FINAL  OPTIMIZATION  INFORMATION 
OBJ  -  . 416732E+00 


114 


39936E+01 


DECISION  VARIABLES  (X-VECTOR) 

1 )  .1 0350E-0 1  . 48287E+00 

CONSTRAINT  VALUES  <G- VECTOR) 

l)  - . 90452E+00  — . 95484E-01  -.67406E-05 

THERE  ARE  1  ACTIVE  CONSTRAINTS 
CONSTRAINT  NUMBERS  ARE 

3 

THERE  ARE  0  VIOLATED  CONSTRAINTS 

THERE  ARE  0  ACTIVE  SIDE  CONSTRAINTS 

TERMINATION  CRITERION 

ABS (OBJ  < I ) —OBJ ( I— 1 ) )  LESS  THAN  DABFUN  FOR  3  ITERATIONS 

NUMBER  OF  ITERATIONS  *  10 

OBJECTIVE  FUNCTION  WAS  EVALUATED  28  TIMES 

CONSTRAINT  FUNCTIONS  WERE  EVALUATED  28  TIMES 

GRADIENT  OF  OBJECTIVE  WAS  CALCULATED  8  TIMES 

GRADIENTS  OF  CONSTRAINTS  WERE  CALCULATED  8  TIMES 


*  C  0  N  M  I  N  * 

*  * 

*  FORTRAN  PROGRAM  FOR  * 

*  * 

*  CONSTRAINED  FUNCTION  MINIMIZATION  * 

*  * 
*************************** 


CONSTRAINED  FUNCTION  MINIMIZATION 
CONTROL  PARAMETERS 


NDV 

ITMAX 

NCON 

NS  IDE 

ICNDIR 

NSCAL 

3 

40 

3 

6 

4 

0 

ITRM 

N1 

N2 

N3 

N4 

NS 

3 

3 

a 

51 

51 

102 

CT 

lOOOOE+OO 

THETA 
. 30000E+0 1 

FDCH 

. lOOOOE— 01 


CTMIN 

.40000E-02 

PHI 

.30000E+01 

FDCHM 
. 10000E-01 


CTL 

lOOOOE-Ol 

DELFUN 
. lOOOOE— 03 

ALPHA X 
. 10000E+00 


LOWER  BOUNDS  ON  DECISION  VARIABLES  (VLB) 

1)  . lOOOOE— 05  . 10000E-05  . tOOOOE-OS 

UPPER  BOUNDS  ON  DECISION  VARIABLES  (VUB) 

l)  . tOOOOE+11  . 10000E+ I I  .  10000E+U 

ALL  CONSTRAINTS  ARE  NON-LINEAR 


CTLMIN 
. 10000E-02 

DABFUN 
. 28184E— 03 

AB0BJ1 
. 10000E+00 


INITIAL  FUNCTION  INFORMATION 

OBJ  -  . 281844E+00 

DECISION  VARIABLES  (X- VECTOR) 

1)  . 12750E-01  • 26953E+00  .39279E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  -.59846E+00  — . 40154E+00  -.89384E-01 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  >  1  OBJ  *  . 2808SE+00 

DECISION  VARIABLES  <X -VECTOR) 

1)  .12705E-01  . 26953E+00  .39279E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

l)  — . 62202E+00  — • 37798E+00  -.90431E-01 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 

ITER  -  2  OBJ  -  . 27999E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 12666E-01  .26953E+00  .39279E+01 

CONSTRAINT  VALUES  (G- VECTOR) 

1)  — . 64234E+00  -.35766E+00  -.91158E-01 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  »  3  OBJ  -  . 277&4E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 12560E-01  . 26932E+00  .39279E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  — • 69737E+00  -.30243E+00  -.80412E-01 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  -  4  OBJ  -  .27733E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  • I2547E— 01  . 26952E+00  .39279E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

I)  — . 70438E+00  -.29562E+00  -.55992E-01 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  «  3  OBJ  -  .27719E+00 

DECISION  VARIABLES  (X-VECTOR) 

t)  . 12340E-01  . 26952E+00  .39279E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  -. 70820E+00  -.29180E+00  -.11972E-02 

*4  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER 


6 


OBJ 


27719E+00 


DECISION  VARIABLES  (X-VECTOR) 

1>  . 12340E-01  .26932E+00  .39279E+01 

CONSTRAINT  VALUES  (Q- VECTOR) 

1)  --70827E+00  - . 29 1 73E+00  -.11783E-04 

**  CONSTRAINT  3  HAS  ZERO  6RADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  -  7  OBJ  -  .27719E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 12540E— 01  . 26932E+00  .39279E+01 

CONSTRAINT  VALUES  <G- VECTOR) 

1)  -.70827E+00  -.29173E+00  16167E-06 


FINAL  OPTIMIZATION  INFORMATION 
OBJ  »  . 277 187E+00 


DECISION  VARIABLES  (X-VECTOR) 

1)  . 12540E— 01  . 26952E+00  .39279E+01 

CONSTRAINT  VALUES  (G- VECTOR) 

1)  -.70827E+00  — . 29 1 73E+00  16167E-06 

THERE  ARE  1  ACTIVE  CONSTRAINTS 
CONSTRAINT  NUMBERS  ARE 
3 

THERE  ARE  O  VIOLATED  CONSTRAINTS 
THERE  ARE  0  ACTIVE  SIDE  CONSTRAINTS 
TERMINATION  CRITERION 

ABS(OBJ ( I ) -OBJ ( 1-1 ) )  LESS  THAN  DABFUN  FOR  3  ITERATIONS 
NUMBER  OF  ITERATIONS  «  7 

OBJECTIVE  FUNCTION  WAS  EVALUATED  22  TIMES 

CONSTRAINT  FUNCTIONS  WERE  EVALUATED  22  TIMES 

GRADIENT  OF  OBJECTIVE  WAS  CALCULATED  7  TIMES 

GRADIENTS  OF  CONSTRAINTS  WERE  CALCULATED  7  TIMES 


C  Q  N  n  I  N 


FORTRAN  PROGRAM  FOR 
CONSTRAINED  FUNCTION  MINIMIZATION 

************************* 


CONSTRAINED  FUNCTION  MINIMIZATION 


CONTROL  PARAMETERS 


IPRINT 

2 

NDV 

3 

ITMAX  NCON 

40  3 

NS IDE  ICNDIR 

6  4 

NSCAL  NFDG 

0  1 

LINOBJ 

0 

ITRM 

3 

N1  N2 

3  8 

N3  N4 

51  51 

N5 

102 

CT 

-. 10000E+01 

CTMIN 
. 40000E— 02 

CTL 

— • IOOOOE— 01 

CTLMIN 
. IOOOOE— 02 

THETA 
. 30000E+0 1 

PHI 

.30000E-M51 

DELFUN 
•  1 OOOOE— 03 

DABFUN 
. 18345E-03 

FDCH 

. IOOOOE— 01 

FDCHM 
• IOOOOE— 01 

ALPHA X 
.  10000E-MX) 

AB0BJ1 
. 10000E+00 

LOWER  BOUNDS  ON  DECISION  VARIABLES  (VLB) 

1)  . IOOOOE— 05  . 10000E-05  .10000E-05 

UPPER  BOUNDS  ON  DECISION  VARIABLES  (VUB> 

1)  . 10000E+1 1  . 10000E+1 1  . 10000E+1 1 

ALL  CONSTRAINTS  ARE  NON-LINEAR 


INITIAL  FUNCTION  INFORMATION 

OBJ  -  . 1 85449E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 1S300E-01  . 15000E+00  .38700E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  50768E+00  -.49232E+00 


99915E-01 


ITER 


1 


OBJ 


1710SE+00 


DECISION  VARIABLES  (X-VECTOR) 

1)  .  15A47E-01  .  136S0E+00  . 38361E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  — .49714E+00  -.50286E+00  -.55735E-01 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  *  2  OBJ  -  .  17108E-MX>  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 15647E-01  . 13650E+00  . 38361E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  — . 497 1 4E+00  -.50286E+00  -.35735E-01 


ITER  *  3  OBJ  -  . 17I08E+00  NO  CHANGE  IN  OBJ 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 13647E-01  . 13650E+00  .383A1E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  -.49714E+00  - • 50286E+00  -.53733E-01 

**  CONSTRAINT  3  HAS  ZERO  GRADIENT 
DELETED  FROM  ACTIVE  SET 


ITER  -  4  OBJ  »  . 17102E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  . 15642E-01  . 13650E+00  .383A1E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

1)  — . 49962E+00  -.30038E+00  . 10340E-02 


FINAL  OPTIMIZATION  INFORMATION 

OBJ  »  ■ 171013E+00 

DECISION  VARIABLES  (X-VECTOR) 

1)  .  15642E-01  .  13650E+00  .38361E+01 

CONSTRAINT  VALUES  (G-VECTOR) 

l)  -.49962E+00  -.50038E+00  . 10340E-02 

THERE  ARE  1  ACTIVE  CONSTRAINTS 
CONSTRAINT  NUMBERS  ARE 


THERE  ARE  0  VIOLATED  CONSTRAINTS 

THERE  ARE  0  ACTIVE  SIDE  CONSTRAINTS 

TERMINATION  CRITERION 

ABS ( OB J (I ) -OBJ < I - 1 ) )  LESS  THAN  DABFUN  FOR  31 TERAT I ONS 

NUMBER  OF  ITERATIONS  *  4 

OBJECTIVE  FUNCTION  WAS  EVALUATED  10  TIMES 

CONSTRAINT  FUNCTIONS  WERE  EVALUATED  10  TIMES 

GRADIENT  OF  OBJECTIVE  WAS  CALCULATED  3  TIMES 

GRADIENTS  OF  CONSTRAINTS  WERE  CALCULATED  3  TIMES 

THIS  DATA  REPRESENTS  A  REDUCTION  IN  THE  MASS  OF  60.686* 


THESE  ARE  THE  EIGENVALUES  FOR  THE  RECONFIGURED  MODEL 


EIGENVALUE  1  IS: 

89637B97E-02 
-  .  4026337 4E+0  1 


AND  ITS  RIGHT  EIGENVECTOR  IS: 


. 20283846E— I 2 
. 20287406E-03 
. 73890674E-04 
- . 20287408E-03 
7389068 IE-04 
17978023E— 10 
-. 1 OOOOOOOE+O 1 
35469399E+00 
. 1 OOOOOOOE+O 1 
. 35469399E+00 


-.20713803E-12 
12120727E-01 
— . 42990834E-02 
. 12120727E-01 
. 42990834E-02 
— . 1396387 IE- 10 
. 1 7264924E— 08 
— . 32444053E-03 
.OOOOOOOOE+OO 
. 324441 13E-03 


ITS  LEFT  EIGENVECTOR  IS  : 


- • 99686638E— O 1 
. 1 OOOOOOOE+O 1 
.31 033788E+00 
. lOOOOOOOE+Ol 
. 3 1 033788E+00 
39333 199E-01 
. 63 1 73826E-0 1 
639063 18E-03 
. 63 1 75826E-0 1 
639063 18E-03 

EIGENVALUE  2  IS: 

- . 89657897E— 02 
. 4026357 4E+0 1 


. 1 49223 14E+00 
• 24836869E- 1 2 
. 88445692E— 0 1 
— . 24590960E- 1 4 
. 88443692E-0 1 
15906706E— 01 
. 20093084E-0 1 
- . 3679 1595E-04 
. 20093084E— 0 1 
— . 3679 1 593E-04 


AND  ITS  RIGHT  EIGENVECTOR  IS: 


. 20283846E- 12 
. 20287406E-03 
. 7389067 4E-04 
20287 408E-03 
7389068 IE-04 
1 7978025E- 1 0 
1 OOOOOOOE+O 1 
-.33469399E+00 
. 1 OOOOOOOE+O 1 
• 35469399E+00 


.2071 3803E- 1 2 
• 12120727E-01 
. 42990834E-02 
12120727E-01 
42990834E-02 
. 1396387 IE- 10 
-. 17264924E-08 
. 32444033E-03 
. OOOOOOOOE+OO 
324441 13E-03 


ITS  LEFT  EIGENVECTOR  IS  : 


-99686658E-01 
,  lOOOOOOOE+Ol 
,  31033788E+00 
,  lOOOOOOOE+O 1 
,  31 033788E+00 
.  39535 1 99E-0 1 
,  63 175826E-0 1 
,  639063 18E-03 
>  63 1 75826E-0 1 
639063 18E-03 


-.1 49223 14E+00 
— . 24836869E- 1 2 
-.884-  .692E-01 
. 24590960E— 1 4 
— . 8S445692E-0 i 
.15906706E-01 
-.20093084E-01 
• 3679 1595E-04 
-.200930B4E-01 
. 3679 1 595E-04 


EIGENVALUE  3  I Si 

-. 1520 189 1E+0O 
.46979208E+01 


AND  ITS  RIGHT  EIGENVECTOR  IS: 


- & 33983387E— 1 1 
. 688068S7E— 02 
. 5202528 IE-03 
- . 68806857E-02 
5202328 IE-03 
- . 50343469E— 1 1 
-. lOOOOOOOE+Ol 
. 17877366E-01 
. 1 OOOOOOOE+O 1 
17877366E-01 


. 1 1201286E— 1 1 
. 21263747E+00 
-.382221 31 E-02 
— . 21263747E+00 
.382221 3 IE-02 
15770762E-10 
-. 3735721 9E-11 
.3025135 l E-02 
— . 29838035E- 1 4 
-.3025 1551 E-02 


ITS  LEFT  EIGENVECTOR  IS  l 


-. 25502832E-02 
. 1 OOOOOOOE+O 1 
. 335 1 1027E+00 
-. 94879 122E+00 
-. 31910987E+00 
-. 1 3559820E-02 
.  44226740E— 0 1 
— . 82829267E-03 
— . 40681695E— 01 
. 78781 128E-03 

EIGENVALUE  4  ISt 

-. 13201 89 1E+00 
— . 46979208E+01 


— . 62046723E-02 
. 164621 16E-14 
— . 76983047E— 0 1 
. 553 1 8496E— 02 
. 72996 150E-01 
. 38728660E— 03 
-.1 029969 IE-01 
-. 6282020 IE-04 
. 93731696E-02 
. 62332975E-04 


AND  ITS  RIGHT  EIGENVECTOR  ISt 


33983387E- 1 1 
68806857E-02 


U201286E— 11 
21263747E+00 


.5202328 IE-03  . 3822213 IE-02 

- . 68806857E-02  . 2 1 263747E+00 

- . 5202528 1 E— 03  - . 38222 1 3 1 E-02 

- . 50343469E- 11  .1 5770762E- 1 0 

- . 1 0000000E+0 1  . 373572 1 9E- 1 1 

. 1 7877366E-0 1  - . 3025 1 55 1 E-02 

. 1 0000000E+0 1  . 29838035E— 1 4 

1 7B77366E— 0 1  . 30251551E-02 

ITS  LEFT  EIGENVECTOR  IS  s 
- . 25502852E-02  . 62046723E-02 

. 1 OOOOOOOE+O 1  - . 1 6462 1 1 6E- 1 4 

. 335 1 1 027E+00  . 76983047E-0 1 

94879 122E+00  55318496E-02 

-.3191 0987E+00  - . 72996 1 50E-0 1 

- . 1 3559820E-02  - . 38728660E-03 

. 442267 40E-0 1  . 1 029969 1 E-0 1 

- . 8282926 7E-03  . 6282020 1 E-04 

- . 4068 1 695E-0 1  - . 9373 1 696E-02 

. 7878 1 1 28E-03  - . 62332973E-04 

EIGENVALUE  5  IS: 

-. 25908072E+00 
. 40835030E+01 


AND  ITS  RIGHT  EIGENVECTOR  ISt 

- . 55304663E— 03  - . 2483622 1 E+00 

- . 80832378E— 04  . 1 1 4834 1 OE+OO 

. 2967 0804E— 05  . 4609357 1 E-02 

.1271 2803E-03  . 1 1 589202E+00 

. 1 878485 1 E-05  . 4603 1 533E-02 

. 10000000E+01  17816479E-14 

-. 4624429 1E+00  -. 13552 177E-02 

- . 1 8558945E— 0 1  - . 29379992E-04 

- . 46662383E+00  - . 52720058E-03 

- . 1 8533957E-0 l  - . 33707454E-04 

ITS  LEFT  EIGENVECTOR  IS  * 

- . 76639446E-0 1  - . 1 67359 1 5E+00 

. 1 OOOOOOOE+O 1  -.281 47774E- 1 2 

. 3 1 654856E+00  - . 94705402E-0 1 

. 10000000E+01  — . 30266423E- 1 4 

.31 654856E+00  -.94705402E-01 

- . 428084 1 3E-0 1  . 1 2489566E-0 1 

. 63476509E-0 1  - . 1 9552827E-0 1 

-.671 37883E-03  - . 82860496E-05 

. 65476509E-0 l  - . 1 9552827E-0 1 
-.671 37883E-03  - . 82860496E-05 


EIGENVALUE  6  IS 


25908072E+00 
- . 40835030E+0 1 


AND  ITS  RIGHT  EIGENVECTOR  I St 


-.55304663E-03 
— . 80832378E— 04 
*  29670804E-05 
. 127 1 2803E-03 
. 1 878485 IE-05 
. 1 OOOOOOOE+O 1 
-. 4624429 1E+00 
1 8558943E— O 1 
46662383E+00 
-. 1 8533957E— O 1 


. 2483622 1E+00 
— . 1 1485410E+00 
-. 4609357 IE-02 
1 1589202E+00 
- . 4603 1 533E-02 
. 17S16479E— 14 
. 1 3552 1 77E-02 
. 29379992E— 04 
. 52720058E-03 
. 33707454E— 04 


ITS  LEFT  EIGENVECTOR  IS  s 


-.76659446E-01 
. 1 OOOOOOOE+O 1 
. 3 1 654856E+00 
. 1 OOOOOOOE+O 1 
.31 654856E+00 
428084 13E-01 
. 65476509E— Ol 
-.671 37883E— 03 
. 65476509E— O 1 
- . 67 1 37883E— 03 

EIGENVALUE  7  I St 

-. 1 2842362E+0 1 
. 8307 161 2E+02 


. 1 67359 15E+00 
.28147774E— 12 
. 94705402E-01 
. 30266423E— 14 
. 947 05402E— O 1 
12489566E-01 
. 1 9552827E— 0 1 
. 82860496E— 05 
. 1 9552827E— O 1 
. 82860496E-05 


AND  ITS  RIGHT  EIGENVECTOR  ISt 


58702 105E-05 
-. 1 8605265E-03 
— . 70882689E-04 
-. 1 S643520E— 03 
—.7101 7753E-04 
-. 1 7046570E-0 1 
. 1 OOOOOOOE+O 1 
. 35047023E+00 
. 999839 63E+00 
. 35041336E+00 


. 20529406E-03 
-. 12034930E-0 1 
-. 42177970E-02 
-. 1 2032994E-0 1 
-.42171 103E-02 
- . 75 1 29393E-03 
-. 82216874E— 15 
— . 47 1691 58E-03 
— . 34265667E-04 
— . 48379340E-03 


ITS  LEFT  EIGENVECTOR  IS  * 


1 56586 75E-01 
. 1 OOOOOOOE+O 1 
. 35089629E+00 
. i OOOOOOOE+O 1 
. 35089629E+00 
1 1 80477 6E-04 
. 19973735E-03 
.65 13025 IE-04 
.  19973735E-03 
.65 13025 IE-04 

EIGENVALUE  8  IS: 

1 2842362E+0 1 
- . 8307 161 2E+02 


-. 1 7234424E— 05 
-.75480239E-10 
-.9841 7340E— 06 
. 1 53 1 4432E— 1 4 
— . 984 1 4663E-06 
. 1 8835S37E-03 
-. 1 2032557E— 0 1 
422221 12E-02 
-. 1 2032557E— 0 1 
-. 422221 12E-02 


AND  ITS  RIGHT  EIGENVECTOR  IS: 


-. 58702 105E-05 
-. 1 8605265E-03 
-. 7 0882689E— 04 
-. 1 8643520E— 03 
—.7101 7753E— 04 
-. 17046570E— 01 
. 1 OOOOOOOE+O 1 
. 35047 023E+00 
. 99983963E+00 
. 3504 1 336E+00 


— . 20529406E— 03 
. 1 2034930E— O 1 
. 42177970E— 02 
. 1 2032994E-0 1 
. 42171 103E-02 
. 75 1 29393E-03 
.82216S74E— 15 
. 47169158E— 03 
. 3426566 7E-04 
. 48379340E— 03 


ITS  LEFT  EIGENVECTOR  IS  : 


-. 1 5658675E-0 1 
. 1 OOOOOOOE+O 1 
. 35089629E+00 
. 1 OOOOOOOE+O 1 
.  35089629E+00 
-. 11804776E-04 
. 1 9973735E-03 
.651 3025 IE-04 
. 1 9973735E-03 
.651 3025 IE-04 

EIGENVALUE  9  IS: 

1 3805354E+0 1 
-. 82480 192E+02 


. 1 7234424E-05 
. 75480239E— 10 
. 98417340E-06 
-. 1 53 1 4432E- 1 4 
. 98414663E-06 
-. 1 8835837E-03 
. 1 2032557E-0 1 
. 422221 12E-02 
. 1 2032557E-0 1 
. 422221 12E-02 


AND  ITS  RIGHT  EIGENVECTOR  IS: 


- . 1 54747S6E-0 1  . 24390597E+00 

- . 3392 1 607E-02  -.11 5986 1 5E+00 

. 10534 196E-03  4613741  IE-02 

- . 3487 1 759E-02  -.11711 92BE+00 

. 1 1074735E-03  -. 46065659E-02 

. 1 OOOOOOOE+O 1  - . 1 S226307E- 1 4 

- . 47275093E+00  . 4390 1 672E-0 1 

- . 1 88675 1 8E-0 1  . 765 1671 6E-03 

- . 47735349E+00  . 4458324 1 E-0 1 

- . 1 88396 1 8E-0 1  .741 23524E-03 

ITS  LEFT  EIGENVECTOR  IS  s 

- . 303983 1 6E-06  . 2793B072E-06 

- . 999834 1 OE+OO  -.3611 372 1 E-04 
- . 35463975E+00  -.13311 750E-04 

. 1 OOOOOOOE+O l  . OOOOOOOOE+OO 

. 35469847E+00  . 63802 1 1 5E-06 

- . 45390 1 83E-08  - . 34669883E-08 

- . 20807377E— 03  -.121 1 7204E-0 1 

-.71 675598E-04  - . 42979627E-02 

. 20854 1 82E— 03  . 121 19207E-01 

. 7 1 838934E-04  . 4298672 IE-02 

EIGENVALUE  10  ISs 

-. 1 3805354E+0 1 
. 82480 192E+02 


AND  ITS  RIGHT  EIGENVECTOR  IS: 

- . 1 5474786E— O 1  - . 24390597E+00 

- . 3392 1 607E—02  . 1 15986 15E+00 

. 1 0534 1 96E-03  . 4613741  IE-02 

- . 3487 1 759E-02  .11711 928E+00 

. 1 107 4735E— 03  . 46065659E-02 

. 1 OOOOOOOE+O l  . 1 8226307E- 1 4 

- . 47275093E+00  -.43901 672E-0 1 

-. 188675 18E-01  -. 76516716E-03 

— . 47735349E+00  -. 44583241 E-0 1 

- .  1 88396 1 8E-0 1  -.741 23524E-03 

ITS  LEFT  EIGENVECTOR  IS  i 

-. 303983 16E-06  -. 27938072E-06 

- . 999834 1 OE+OO  .3611 372 1 E-04 

-. 35463975E+00  . 1331 1750E-04 

. 1 OOOOOOOE+O 1  . OOOOOOOOE+OO 

. 3546984 7E+00  -. 638021 15E-06 

- . 45390 1 83E-08  . 34669883E-08 

- . 20807377E— 03  .1211 7204E-0 1 


71675598E-04 
20854 182E-03 
7 1838934E— 04 


. 42979627E— 02 
121 19207E-01 
4298672 IE-02 
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